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ABSTRACT 

We present the galaxy stellar mass function and its evolution in clusters from z ~ 0.8 to the 
current epoch, based on the Wide-field Nearby Galaxy-cluster Survey (WINGS) (0.04 ^ z ^ 
0.07), and the ESO Distant Cluster Survey (EDisCS) (0.4 ^ z < 0.8). We investigate the total 
mass function and find it evolves noticeably with redshift. The shape at M* > 1O U M0 does 
not evolve, but below M* ~ 10 10 *Mq the mass function at high redshift is flat, while in the 
Local Universe it flattens out at lower masses. The population of M» = 10 10 2 - 10 ia8 M Q galax- 
ies must have grown significantly between z = 0.8 and z = 0. We analyze the mass functions 
of different morphological types (ellipticals, SOs and late-types), and find that also each of 
them evolves with redshift. All types have proportionally more massive galaxies at high- than 
at low-z, and the strongest evolution occurs among SO galaxies. Examining the morphology- 
mass relation (the way the proportion of galaxies of different morphological types changes 
with galaxy mass), we find it strongly depends on redshift. At both redshifts, ~ 40% of the 
stellar mass is in elliptical galaxies. Another ~ 43% of the mass is in SO galaxies in local 
clusters, while it is in spirals in distant clusters. To explain the observed trends, we discuss the 
importance of those mechanisms that could shape the mass function. We conclude that mass 
growth due to star formation plays a crucial role in driving the evolution. It has to be accompa- 
nied by infall of galaxies onto clusters, and the mass distribution of infalling galaxies might be 
different from that of cluster galaxies. However, comparing with high-z field samples, we do 
not find conclusive evidence for such an environmental mass segregation. Our results suggest 
that star formation and infall change directly the mass function of late-type galaxies in clusters 
and, indirectly, that of early-type galaxies through subsequent morphological transformations. 

Key words: galaxies: clusters: general — galaxies: evolution — galaxies: formation — galax- 
ies: mass function — galaxies: ellipticals and lenticulars, cD 
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1 INTRODUCTION 

The distribution of galaxy stellar masses at the present day and in 
the past is of fundamental importance for studying the assembly of 
galaxies over cosmic time. 

It is well known that galaxy properties depend both on 
galaxy mass and on the environment the galaxies are part of 
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(e.g. see iBaldrv et alj [20061 : Ivan der Well 120081: IBolzonella et al I 



20091: iTasca et alj|2009h . Several studies (e. g. iBundv et alj|2005l : 



Verea nTet al.ll2008h have shown that mass plays a crucial role in 



determining galaxy properties and in driving their evolution. Color, 
specific star formation rate and in ternal structure are stro ngly cor- 
related with galaxy stellar mass ( Kauffmann et al , 2003-). On the 
other hand, the stellar mass function of galaxies depends on lo- 
cal g alaxy density too (e.g. IBaldrv etal 120061 : IBolzonella et al.l 
2009), showing that the mass distribution is itself related to envi- 
ronment. In addition, the efficiency of environmental mechanisms 
can depend on the mass of the galaxy on which they act, and this 
could cause enviro nmental effects to be interpreted as mass-effects 
dTasca et al.ll2009h . 

Clearly, galaxy mass and environment are strictly linked and 
it is important to analyze them distinctly, trying to separate their 
relative roles. To disentangle mass and environment it is useful to 
study the dependence of some properties on mass fixing the envi- 
ronment, or analyze the importance of the environment fixing the 
mass. In the nature vs nurture scenario, mass represents the pri- 
mary "intrinsic property" closely related to primordial conditions, 
while the environment includes all the possible external processes 
that can influence galaxy evolution. 

Stellar mass assembly is usually characterized by the galaxy 
stellar mass function, that describes the distribution of galaxy 
masses, at different epochs and different environments. Several 
studies have been carried out focusing on th e importance of the 
galaxy mass functions in the general field (e.g . iFontana et al 
Drorv et al.l 120051: iGwvn & Hartwick! 120051 : IBundv et all 



2006 



2006 



Pozzetti et aU2007l) . 



For galaxies with M ^ 10 Mp, several studies (eg. 
IFontana et al] |2004 IBundv et all 120061 : iBorch et all l2006h have 
found that, overall, the evolution of the total mass function from 
z = 1 to z = is relatively modest, which implies that the evolution 
of objects with mass close to the l ocal characteri s tic ma ss is essen- 
tially complete by z~l. However. IFontana et alj Ii2006) found that 
less massive galaxie s evolve more than massive ones. Moreover, 
IPozzetti et al. (2009) showed a very fast rise of the total mass func- 
tion from z~l for M ^ 1O U M0 and they found this trend is due 
to star formation driving the mass assembly at intermediate and 
low masses. Instead, at high masses t his process does not p lay an 
important role and, as a consequence, IPozzetti et ail (2009) found 
again that at high masses there is an almost negligible evolution 
from z ~ 1 to z ~ 0. 

It is interesting to investigate how galaxies of different types 
give a different contribution to the mass function, shaping, for ex- 
ample, either the massive tail or the low mass tail of the total mass 
distribution. There are several ways to subdivide galaxies into at 
least two populations, usually either according to their star for- 
mation histories (being passive or star-forming, for example us- 
ing a rest-frame color, hence separating blue and red galaxies, their 
SEDs, their spectroscopic features) or their structure (structural pa- 
rameters or morphologies). 

In the Local Universe, a sort of bimodality in galaxy proper- 
ties has been observed: more massive galaxies tend to be passive, 
have high surface mass densities, redder colors and have spectra 
characterized by a pronounced break at A -~4000 A, while lower- 
mass galaxies have low surface mass densities, tend to be still 
star forming, sh ow bluer colors and have spectra characterized b y 

120031 ~ 



emission lines (Kauffmann et al. 



Baldr v et alj d2004l) . IBaldrv et al 



.Bri nchmann et al. 2004). 
(20061) and IBaldrv et alj 12008). 



subdividing galaxies depending on their colors, found a bimodal 
shape in the local mass function in the field with an upturn related 



to the two different populations: early-type galaxies dominate the 
high masses, while l ate-types dominate at low/intermediate masses 
dPozzetti et alj2009h . 

Many studies have been performed to investigate the be- 
havior of t his bimodality at higher redshifts: for example, 
IBorch et al. (200^), who s ubdivide galaxies d epen ding on their 
color IFontana et alj J2004h . IBundv et alj J2006l) and lVergani et al.1 
(2008), who use spectral features to subdivide galaxies with differ- 
ent properties. Although using different techniques, all of them em- 
phasize that the mass functions separately for early-type and late- 
type galaxies depend on redshift. Early-type galaxies always dom- 
inate the massive end, while late-type galaxies mostly contribute to 
the intermediate/low-mass part of the mass function at all redshifts. 
From z ~ 1 to z ~ 0, early-type galaxies show a rapid evolution in 
their mass function, while late-type galaxies show a relatively little 
change. This has been interpreted as the result of mass growth due 
to star formation in late-type galaxies, accompanyed by a subse- 
quent cessation of star formation in some of the blue galaxies, that 
turn red and con tribute to the build-up of the red galaxy population 
dBell et al.ll2007l) . 

Few studies have attempted to investigate the mass functions 
of different galaxy types as a function of environment ( s ee e 
iBaloghetaimOOll : IBundv et alj l2006h - IBolzonella etal] dioO 
classifying galaxies following a spectral energy distribution clas- 
sification, found that the shape of the mass function of different 
galaxy types is the same in environments with different local den- 
sity, so they argued for the existence of a quasi-universal mass func- 
tion for each type regardless of environment, and also argued for a 
quasi-universal evolution of the mass function of each type. 

The studies mentioned above all focus on the mass functions 
of galaxies with different stellar histories in the general field, but a 
clear picture on the stellar mass assembly and how it depends on 
the mass (mass-assembly downsizing) and galaxy types has not yet 
emerged. 

Fewer studies have focused on the mass functions as a 
function of galaxy morphology. All of t hem consider g alaxie s 
in the general fie l d. IBundv et alj d2005l) : IPannella et all d2006h : 
iFranceschini et al] d2006h investigated mass functions of different 
morphological types, using ACS/HST images to classify galaxies. 
They found that, at z ~1, morphologically early- and late-types 
are present in similar numbers at all masses. By z ~ 0.3, ellipti- 
cals dominate the high mass population, suggesting that merging 
or some other transformation process is active. At all redshifts, they 
found that late-type galaxies dominate at lower masses, while early- 
type galaxies become prominent at higher masses. 

Using both a spectral and a morphological classification, 
IllbertetalJ d2010h found that most of massive quiescent galaxies 
have an elliptical morphology at z <0.8. They also show similar 
mass functions. Therefore, a dominant mechanism has to link the 
shutdown of star formation and the acquisition of an elliptical mor- 
phology in massive galaxies. 

As the evolution of the mass functions of passive versus star- 
forming galaxies must be influenced by the quenching of star for- 
mation in previously star-forming galaxies (Bell et al. 2007), the 
evolution in the mass functions of different morphological types 
must be influenced by the transformation between z = 1 and z = 
of some late-type galaxies into early-type galaxies, of which evi- 
dence has accum ulated both in the fi eld and, especially, in clusters. 

In the field, IBundv et all d2005l) found evidence for the trans- 
formation of late- into early-type galaxies as a function of time. 
Based on the observed mass functions, this transformation process 
appears to be more important at lower masses (M, ^ 1O H M0 ) 
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because the most massive early - type g alaxies are already in place 
at z~l. Similarly. lOesch et alj ^20 id) showed that the evolution 
of the morphological distribution in galaxies since z ~ 1 depends 
strongly on stellar mass: ellipticals dominate the entire population 
since z=l, and the Hubble sequence has quantitatively settled down 
by this epoch. Most of the morphological evolution from z=l to 
z=0.2 takes place at masses M» ^ 10 u Mq, where the fraction of 
spirals substantially drops and the contribution of early-types in- 
creases. 

ICapak etail d2007h found that in the field galaxies are trans- 
formed from late- to early-type galaxies more rapidly in dense than 
in sparse regions. Comparing their field data to cluster data, they 
found that in clusters transformations occur even more rapidly. 

The morphological evolution is much more striking in clus- 
ters, where indeed was discovered: the fraction of late-type galax- 
ies is much higher, and the fraction of SO galaxies proportionally 
lower, in distant than in nearby clusters, suggesting an evolution 
from late to early-type galaxies for a significant fraction of today's 
early-type galaxies in clusters (Dressier et al. 1997, Fasano et al. 
2000, Postman et al. 2005, Smith et al. 2005, Desai et al. 2007). 
Interestingly, the same studies find that the fraction of ellipticals in 
clusters does not evolve with redshift. 

It isn't surprising that the impact of morphological transfor- 
mations depends on environment. In clusters, galaxies are expected 
to be affected by a number of physical processes that can al- 
ter the course of their evolution. Galaxies in dense environments 
are subject to a variety of external stresses, which are, in gen- 
eral, not cond uctive to the maintenance of spi ral structure: i.e. 
ram pressure dGunn & Gottl 1 19721 : iBekkil 12009^ galaxy harass- 
ment dMooreetal.lll996l) . cluster tidal forces dBvrd & Valtoneq 
1990), interaction/merging ( Icke 1985; Bekki 19 9jl), stran gulation 
dLarson et alJl98CtlFont et alj2008l : lMcCarthv et ajj2008t) can act 
with different efficiency depending on environment, so field galax- 
ies infalling into larger structures can be transformed from gas-rich 
spirals to gas-poor lenticular galaxies. 

Until now, a few studies of the mass function of galaxies in 
clusters has been c a rried o ut, due to the few observations available. 
Kodama & Bower (2003) have estimated the stellar mass func- 
tion for two clusters at z ~ 1 transforming their K s -band lumi- 
nosity function by correcting it for the star formation contribu- 
tion as estimated from the J-K j colours. This and other works (e.g . 



Aragon-Salamanca et al.lll993r.lBarger et alJl998l ; rToft et al-fcool 
Andreonll2006l : lDeProprisetalJl2007h have highlighted the fact 



that the evolution of the characteristic luminosity at red wave- 
lengths in clusters is consistent with passive evolution and with 
massiv e cluster galaxie s being fully assembled by z ~ 1. In par- 
ticular, iToft et al j d2004l) pointed out an evolution in the faint end 
slope of the K,-band luminosity function of one cluster at z = 1.2, 
suggesting that clusters at z ~ 1 contain relatively smaller fractions 
of low mass galaxies than clusters in the local Universe. 

In this context, the Wide -field Nearby Galaxy-cluster Survey 
(WINGS) dFasano et all2006l) and the ESO Distant Cluster Survey 
(EDisCS) (White et al. 2005), two surveys whose main aims are 
to characterize clusters and their galaxies at low and high redshift 
respectively, represent a significant increase in the number of well 
studied clusters at low and high redshift respectively. 

In this paper we focus on the evolution of the galaxy mass 
functions in clusters. We first analyze the total mass function, and 
then separately the mass functions of different morphological types 
and the evolution of the morphological fractions, with the aim of 
understanding the mechanisms that drive the evolution of the mass 
distributions. 



This paper is organized as follows. In §2 we present the data- 
sets used (WINGS and EDisCS), describing the surveys and the 
clusters selection at different redshifts; in §3 we derive stellar 
masses. In §4 we present the mass and magnitude limited galaxy 
samples, and other galaxy properties we need for the analysis. We 
then show the total mass function and its evolution with redshift 
(§5.1), and the mass functions of different galaxy morphological 
types (§5.2) and their evolution from z ~ 0.8 to the current epoch 
(§5.3). In §6 we analyze the morphology-mass relation, the evolu- 
tion of the morphological fractions and of the mass fractions in each 
morphological type. In §7 we comment our findings, discussing the 
possible mechanisms that are responsible for the observed evolu- 
tion. Finally, in §8 we summarize our conclusions. 

Throughout this paper, we assume Ho = 70kms~' Mpc~', O m = 
0.30, Oa = .70. The adopted initial mass function (IMF) is a 
iKroupal d200lh in the mass range 0.1-100 M @ . 



2 THE DATA SAMPLES 

To perform the study of the mass distribution of galaxies of dif- 
ferent morphological types and its evolution from z ~ 0.8 to the 
current epoch, we assemble two galaxy cluster samples in two 
redshift intervals: the sample at low-z is sele cted from the Wlde - 
field Nearby Galaxy-cluster Survey (WINGS) dFasano et alj200d) . 
while the sample a t high-z is selected from the ESO Distant Cluster 
Survey (EDisCS) dWhite et al.ll2005l) . 



2.1 Low-z sample: WINGS 

WINGS3 ( Fasano et al.|[2006l) is a multiwavelength survey of clus- 
ters at 0.04 < z < 0.07 whose main goal is to characterize the pho- 
tometric and spectroscopic properties of galaxies in nearby clusters 
and to describe the changes of these properties depending on galaxy 
mass and environment. 

Clusters were selected in the X-ra y from the ROSAT Bright- 



est Cluster Sample and its extension (Ebel ing et al 



11998112001 

and t he X-ray Brightest Abell-type Cluster sample dEbeling et al 
1996). WINGS clusters cover a wide range of velocity dispersion 
Ccius (typically 500- 1 100 kms~' ) and a wide range of X-ray lumi- 
nosity L x (typically 0.2-5 x 10 44 erg^'). 

The surve y is mainly based on optical B, V imaging of 78 
galaxy clusters dVarela et al.l2009l) . that has been complemented by 
a spectroscopic survey of a subsample of 48 clusters, ob tained with 
the sp ectrographs WYFFOS@WHT and 2dF@AAT dCava et al.1 
2009), by a near-infrared (J, K) sur vey of a subsample of 28 clusters 
obtained with WFC AM @ UKIRT dValentinuzzi et alj2009l) . and by 
U broad-band and H a narrow-band imaging of a subset of WINGS 
clusters, obtained with wide-field cam eras at different telescopes 
(INT, LBT, Bok) dOmizzolo et alj20ld) . 

The spectroscopic target selection was based on the WINGS 
B, V photometry. The aim of the target selection strategy was to 
maximize the chances of observing galaxies at the cluster redshift 
without biasing the cluster sample. We selected galaxies with a total 
V ^ 20 magnitude and a V magnitude within the fiber aperture of V 
< 21.5 and with a color within a 5 kpc aperture of (B-V)skpc ^ 1.4, 
to reject background galaxies, much redder than the cluster red se- 
quence. The exact cut in color was varied slightly from cluster to 
cluster in order to account for the redshift variation and to optimize 
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cluster name z a DM 

(kms~ l ) (mag) 



A 1069 


0.0653 


690± 68 


37 


.34 


A119 


0.0444 


862± 52 


36 


.47 


A151 


0.0532 


760±55 


36 


.87 


A1631a 


0.0461 


640±33 


36 


.55 


A1644 


0.0467 


1080± 54 


36 


.58 


A2382 


0.0641 


888± 54 


37 


.30 


A2399 


0.0578 


712± 41 


37 


.06 


A2415 


0.0575 


696±51 


37 


.05 


A3 128 


0.06 


883± 41 


37 


.15 


A3158 


0.0593 


1086± 48 


37 


.12 


A3266 


0.0593 


1368± 60 


36 


.12 


A3376 


0.0461 


779± 49 


36 


.55 


A3395 


0.05 


790± 42 


36 


.73 


A3490 


0.0688 


694± 52 


37 


.46 


A3556 


0.0479 


558± 37 


36 


.64 


A3560 


0.0489 


710± 41 


36 


.68 


A3809 


0.0627 


563± 40 


37 


.25 


A500 


0.0678 


658±48 


37 


.42 


A754 


0.0547 


1000± 48 


36 


.94 


A957x 


0.0451 


710± 53 


36 


.50 


A970 


0.0591 


764± 47 


37 


.11 



Table 1. List of WINGS clusters analyzed in this paper and their redshift z, 
their velocity dispersion a and their distance modulus DM. 

the observational setup. These very loose selection limits were ap- 
plied so as to avoid any bias in the colors of selected galaxies. 

In this paper, we only consider spectroscopically confirmed 
members of 21 of the 48 clusters with spectroscopy. This is the 
subset of clusters that have a spectroscopic completeness (the ra- 
tio of the number of spectra yielding a redshift to the total number 
of galaxies in the photometric catalog) larger than 50%. The com- 
pleteness is determined using V-band magnitudes and turns out to 
be rather flat with magnitud e for most clusters, for the reasons dis- 
cussed in lCava et alj d2009). Moreover, completeness is essentially 
independent from the distance to the centre of the cluster for most 
clusters. 

Our imaging covers a 34' x 34' field, that corresponds to about 
0.6 i?200, for all the 21 clusters considered, except for A1644 and 
A3266 where the ~ 0.5/?200 is covered. i?200 is defined as the radius 
delimiting a sphere with interior mean density 200 times the critical 
density of the Universe at that redshift, and is commonly used as an 
approximation for the cluster virial radius. The ^200 val ues for our 
structu res are computed from the velocity dispersions bv lCava et al.l 
d2009h . 

Galaxies are considered members of a cluster if their spec- 
troscopic redshift lies within ±3<r from t he cluster mean r edshift, 
where a is the cluster velocity dispersion dCava et alj|2 009). 

The clusters used in this analysis are listed in TableQ] 

2.2 High-z sample: EDisCS 

EDisCS is a multiwavelength photometric and spectroscopic sur- 
vey of galaxies in 20 fields containing galaxy clusters at 0.4 < z < 1 
(White etalj2005h . Its main goal is to characterize both the clusters 
themselves and the galaxies within them. 

Clusters were drawn fr om the Las Campan as Distant Clus- 
ter Survey (LCDCS) catalog dGonzalez et alJuOOll) . They were se- 
lected as surface brightness peaks in smoothed images taken with 
a very wide optical filter (~ 4500-7500 A). The redshifts of the 



LCDCS cluster candidates were initially estimated from the ap- 
parent magnitude of the brightest cluster galaxy (BCG). The 20 
EDisCS fields were chosen among the 30 highest surface bright- 
ness candidates, after confirmation of the presence of an apparent 
cluster and of a possible red seq uence with VLT 20 min exposures 
in two filters (White et alj[2005l) . Then, all candidates were subdi- 
vided in two groups, the first one at intermediate redshift (z ~ 0.5) 
and the second one at high redshift (z ~ 0.8). 

For all 20 fields, EDisCS has obtained deep op tical multi- 
band photometry with FORS2/V LT dWhite et alj|2005h and near - 
IR photometry with SOFI/NTT dAragon-Salamanca et al1l2010h . 
Photometric red shifts were mea sured using both optical and in- 
frared imaging dPello et alj|2009h . Their determinations, their per- 
formance, and their use to iso late cluster members, are described 
in detail in IPello et alj d2009h and iRudnick et ail d2009h . and are 
briefly summarized in the following. Photometric redshifts were 
computed for every object in the EDisCS fields using two inde- 
pende nt codes, a modified ve rsion of the pub licly available Hyperz 
code dBolzonella et alj200ol) and th e code oflRudnick et"ai] ( 2001) 
with the modifications presented in lRudnick et alj d2003h . The ac- 
curacy of both methods is a(Sz) ~ 0.05 -0.06, where Sz = zs \ e +~J e h °' ■ 
Membership was establish ed using a modifi e d vers ion of the tech- 
nique first developed in iBrunner & Lubinl (2000), in which the 
probability of a galaxy to be at redshift z (P(z)) is integrated in a 
slice around the cluster redshift to give P c iust for the two codes. The 
width of the slice around which P{z) is integrated should be of the 
order of the uncertainty in redshift for the galaxies in question. In 
our case a Az = ±0. 1 slice around the spectroscopic redshift of the 
cluster Zdust was used. A galaxy was rejected from the member- 
ship list if P c i ust was smaller than a certain probability P t i, res h for 
either code. The Pthresh value for each cluster was calibrated from 
our spectroscopic redshifts and was chosen to maximize the effi- 
ciency with which we can reject spectroscopic non-members while 
retaining at least ~ 90% of the confirmed cluster members, inde- 
pendent of their rest-frame (B-V) color or observed (V-I) color. In 
practice we were able to choose thresholds such that we satisfied 
this criterion while rejecting 45%-70% of spectroscopically con- 
firmed non-members. Applied to the entire magnitude limited sam- 
ple, our thresholds reject 75%-93% of all galaxies with I tot < 24.9. 
A posteriori, we verified that in our sample of galaxies with spec- 
troscopic redshift and above the mass limit described below, ~20% 
of those galaxies that are photo-z cluster members are spectroscop- 
ically interlopers and, conversely, only ~6% of those galaxies that 
are spectroscopic cluster members are rejected by the photo-z tech- 
nique. 

Deep spectroscopy with FORS2/VLT was obtained fo r 18 of 
the fields dHallidav et alJl2004l ; lMiivang- Jensen et alj|2008l) . Spec- 
troscopic targets were selected with the aim to produce an unbiased 
sample of cluster galaxies. Conservative rejection criteria based 
on photometric redshifts (different from those adopted to assess 
the photo-z membership described above) were used in the selec- 
tion of spectroscopic targets to reject a significant fraction of non- 
members, while retaining an unbiased spectroscopic sample equiv- 
alent to a purely I-band selected one. Typically, spectra of more 
than 100 galaxies per field were obtained, with 1^22 for the z~0.5 
cluster candidates and 1^23 for the z~ 0.8 ones. 

ACS/HST mosaic imaging in F81 4W of 10 of the h ighest red- 
shift clusters has also been acquired dDesai et alj|2007l) . covering 
with four ACS pointings a 6.5' x 6.5' field with an additional deep 
pointing in the centre. Thi s field covers the R200 of all clusters, ex- 
cept for cl 1232.5-1250 dPoggianti etai]|2006l) . The R 2 oa values 
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cluster name 



(kins ) 



cl 1040 
cl 1054 
cl 1054 
cl 1103 
cl 1138 
cl 1216 
cl 1227 
cl 1232 
cl 1354 



7-1155 

4- 1146 
7-1245 

7- 1245 
2-1133 

8- 1201 

9- 1138 

5- 1250 
2-1230 



0.70 
0.70 
0.75 
0.62 
0.48 
0.79 
0.64 
0.54 
0.76 



418 
589 



504 



-70 
-113 



-65 

336* 

732 +72 



101 



-76 

4 

574^ 
1080+J' 9 
64*1 



Table 2. List of EDisCS clusters analyzed in thi s paper, with cluster 
name, redshift z and vel ocity dispersion a (from Halliday et al] |2004| ; 
Milva ng- Jensen et alj2008T) . 



for our structures are computed from the velocity dispersions by 
IPoggianti etaD J2008h . 

In this paper, to have a large sample of galaxies at high redshift 
(see below), we analyze the data of photo-z members of 9 EDisCS 
clusters, for which HST images and, consequently, morphologies 
were available. Even though HST data for cl 1037.7-1243 are avail- 
able, this cluster has no galaxies above the completeness limits we 
will adopt, hence in our analysis we do not consider it. Table [2] 
presents the list of clusters used. 



3 GALAXY STELLAR MASSES ESTIMATES 



In this paper we adopt the iKroupal j20(5Th IMF in the mass range 
0.1-100 Mq to compute the galaxy stellar masses . We determine 
them i n the same way for both data-sets. We follow Bel l & de Jond 
who used a spectrophotometric model finding a strong cor- 
relation between stellar mass-to-light (M/L) ratio and opticals col- 
ors of the integrated stellar populations for a wide range of star 
formation histories. They found that this correlation is quite robust 
to uncertainties in stellar population and galaxy evolution models. 

We determine stellar masses using the relation between 
M/Lb and rest-frame (6 - V) color and the equation given in 
iBell&de JorielfeOOlh : 

log 10 (M/L B ) = a B + b B (B-V) (1) 

For the Bruzual & Chariot model with a lSalpeteri d 19551) IMF (0. 1- 
125 Mq) and solar m etallicity , a B = - 0.51 and b B = 1.45. Then, we 
scale our masses to a Kroupal d200lh IMF adding -0.19 dex to the 
logarithmic value of the masses. 

The choice of the photometric bands is dictated by the fact 
that we have only B and V for all the WINGS clusters used, and 
also that the rest-frame B-band has been directly observed for all 
EDisCS clusters. K-band luminosity would be preferable to deter- 
mine masses, as the K-band is less sensitive to recent star forma- 
tion and dust. However, K-band is available only for a subset of 
the WINGS clusters considered, and K-band wavelengths have not 
been directly observed by EDisCS, for which K-band magnitudes 
were extrapolated from the photo-z templates. 

3.1 Low-z: WINGS 

For WINGS spectroscopically confirmed cluster galaxies, the total 
luminosity L B has been derived from the t otal (Sextractor AUTO) 
observed B magnitude jVarela et al.1l2009l) . corrected for distance 
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Figure 1. Comparison between the determination of WINGS masses 
used in this paper and those derived from the spectr ophotometric m odel 
<Fritzetalj2 007. 2010). All masses are converted to a Kroupa] <200ll) IMF 
in the mass range 0.1-100 Mq. Points are logarithmic values. The blue line 
is the best fit to the points, and the red line is the 1 : 1 relation. 



modulus and foregrou nd Galaxy extin ction, and k-corrected using 
tabulated values from Poggianti U1997I) . The B-V color used to cal- 
culate masses was derived from observed B and V aperture mag- 
nitudes measured within a diameter of 10 kpc around each galaxy 
baricenter, corrected as the total magnitude. 

Stellar masses for WINGS galaxies had been previously deter- 
mi ned by fitting the o ptical spectrum (in the range ~ 3600- ~ 7000 
A) dFritz et al.l201fj|). with th e spectro-photometric model fully de- 



scribed in Fritz et ail {2007). In this model, all the main spectro- 
photometric features are reproduced by summing the theoretical 
spectra of Simple Stellar Population (SSP) of 13 different ages 
(from 3 x 10 6 to ~ 14 x 10 9 years). Dust extinction is allowed to 
vary as a function of SSP age, and the metallicity can vary among 
three values: Z=0.004, Z=0.02 and Z=0.05. These mass estimates 
were then corrected for color gradients within each galaxy (i.e. for 
the difference in color within the fibre and over a lOkpc diame- 
ter). For a detailed descr iptio n of the de t ermin ation of masses see 
IValentinuzzi et alj ( 1201 Oh and lFritz et al] ( |2010|) . 

The comp arison of masses u sed in this paper with masses de- 
termined by the lFritz et al.1 {2010) spectro-photometric model con- 
verted to the adopted IMF (see Fig.QJ shows no offset between the 
two estimates, and an rms of ~ 0.3 dex that we adopt as the typical 
error on the masses. Above logM»/M = 9.8, that is the WINGS 
mass completeness limit (see 34.1. Q there is an offset of ~ 0.06 
dex and the rms is ~ 0.3 dex. An exhaustive comparison among 
masses and error s determined from different bands can be found in 
iFritz etaljjlolch . 

For homogeneity with the EDisCS m ass estimates (see nex t 
section), throughout this paper we only use iBell & de Jonl feOOlh 
masses. 
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Figure 2. Comparison between the determination of EDis CS masses used 
in this paper and those derived using the k-correct tool iBlanton & Roweisl 
|2007|) . Only galaxies for whi ch spectroscopy is available are considered. 
All masses are converted to a Kroupa (2001) IMF in the mass range 0.1- 
100 Mq. Points are logarithmic values. The blue line is the best fit to the 
points, and the red line is the 1:1 relation. 

3.2 High-z: EDisCS 

For EDisCS galaxies, stellar masses for spectroscopic members 
were estimated by John Moustakas (2009, privat e communication) 
using the kcorrect tool (Blanton & Roweis 2007fl that models the 
available observed broad band photometry, fitting templates ob- 
tained with spectrophotometric models. 

For this paper, to be fully consistent with stellar masses of the 
sample at low-z, we ha ve estimated stellar ma sses following again 
the method proposed by Bell & de Jong ( 2001), using total absolute 
magnitudes derived from photo-z fitting fixing each galaxy redshift 
to be equal to the spectroscopic redshi ft of the cluster it belongs to 
dRudnick etal.| [2003: IPello et al]|2009l) . converting our values to a 
iKroupal fcOOlh IMF. 

In Figure [2] a comparison of the two methods is shown only 
for spectroscopic members of all EDisCS clusters (not just the 
9 with HST used in this paper). There is a satisfactory agree- 
ment, with a systematic offset of 0.03 dex and an rms of 0.08 
dex. If we consider data only above the EDisCS completeness limit 
(logM»/M Q ^ 10.4, see 94.2.11 ), we find a mean offset of ~ 0.09 
dex and an rms of 0.08 dex. 

Finally, we note that a only slightly larger systematic offset 
and rms (0.13 and 0.13, respectively) are found when comparing 

the masses of EDisCS spectroscopic members derived with the 

I 1 I k 

Bell & de Jons (2001) method with those derived b y fitting the 
EDisCS spectra using the non-parametric method of [ Ocvirk et al.l 
d2006h : IOcvirkl feoiolb . t he new models oflVazdeklsl fcq09l) and the 
MILES stellar library of lSanchez-Blazquez et al. (2006). 



2 http://cosmo.nyu . edu/mb 1 44/kcorrec t/ 



4 GALAXY SAMPLES 

Both at high- and low-z, we analyze the mass functions and the 
fractions of galaxies of different morphological types, as deter- 
mined by high-quality imaging. For this purpose, galaxies are di- 
vided in three groups: ellipticals, lenticulars (or SOs) and late-type 
galaxies (spirals + irregulars). In some cases, we will also consider 
ellipticals and SOs together, defined as early-type galaxies. 

In the following, depending on our aim, we use two different 
methods of galaxy sample selection. 

To study the mass distribution of different morphological 
types (§5), the best is to use a mass limited sample. In this case, 
the mass limit endorsed needs to be such to ensure completeness, 
i.e. to include all galaxies more massive than the limit regardless 
of their color or morphological type. To determine this limit, we 
compute the mass of an object whose observed magnitude is equal 
to the faint magnitude limit of the survey, and whose color is the 
reddest color of a galaxy at the highest redshift considered. 

In the subsequent part of the paper (§6), we compare the evo- 
lution of the mass distributions with the evolution of the morpho- 
logical fractions. Since in the literature the latter are available only 
for a magnitude limited sample, we assemble also galaxy magni- 
tude limited samples. In this case, the minimum detected stellar 
mass depends both on redshift and galaxy color. 



4.1 Low-z sample: WINGS 

For WINGS, we consider only spectroscopically confirmed mem- 
bers of our clusters. Using all galaxies in the photometric sample 
would provide a larger dataset, however for WINGS the photomet- 
ric data cannot be used to assess cluster membership due to the low 
redshift and restricted photometric coverage. 

Since galaxies with available spectroscopy are only a fraction 
(> 50% in the clusters used) of the total number of possible spec- 
troscopic targets, we need to apply a statistical correction to correct 
for incompleteness. This is obtained weighting each galaxy by the 
inverse of the ratio of the number of spectra yielding a redshift to 
the tot al number of gala xies in the photometric catalog, in bins of 
1 mag JCava et alj|2009l) . 

As in several other works (see e.g. lvon der Linden et al.l2O10h . 
in each cluster, we exclude the BCG, defined as the most luminous 
galaxy of each cluster, that could alter the mass distribution of the 
whole sample. In fact, the properties of BCGs are in many respect 
peculiar with respect to other galaxies, and they are the subject of 
many se parate studies dedicated only to this class of objects. For 
example. lFasano et al.l (2010) found that also in the WINGS dataset 
BCGs are heterogeneus with respect to the global population of el- 
lipticals for their distribution of apparent (and intrinsic) flattenings. 

Only galaxies lying within O.6.R200 are considered, because 
this is the largest radius covered approximately in all clusters. 

Morphological types are derived from V-band images using 
MORPHOT, an automatic tool for galaxy morphology, purposely 
devised in the framework of the WINGS project. MORPHOT was 
designed with the aim to reproduce as closely as possible visual 
morphological classifications. 

MORPHOT extends the classical CAS (Concentra - 
tion/Asymmetry/clumpinesS) parameter set dConselicd 120031) . 
by using 20 image-based morphological diagnostics. Fourteen of 
them have never been used, while the remaining six [besides the 
CAS parameters, the Sersic index, the Gini and M20 coefficients 
jLotzetalJ 120041) 1 are actually already present in the literature, 
although in slightly different forms. Defining the newly introduced 
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Figure 3. Normalized distributions of the MORPHOT Azimuthal Coeffi- 
cient for the visually classified ellipticals (366 objects, red histogram) and 
SO galaxies (267 objects, blue histogram) of the MORPHOT calibration 
sample. 



o 
t/1 











43% 

o 
o 






16% 


90% 

CO 


57% 




18% 

CN 


56% 

in 


10% 

M 

CN 




36% 

c\ 


58% 

m 


24% 

N 






64% 

to 


24% 

o 
r-O 


4% 







E SO SpE SpL Irr 

Visual Type 



Figure 4. Comparison between visual and MORPHOT broad morphologi- 
cal types for the 93 1 galaxies of the MORPHOT calibration sample. In the 
2D bins of the plot the percentages of visual/MORPHOT broad types in 
different MORPHOT/visual types are reported in the top/right sides of the 
bins. 



diagnostics and explaining their meaning i s beyond the scope 
of the present paper. We refer the reader to iFasano et 
Appendix A therein) for an outlining of the logical sequence 
and the basic procedures of MORPHOT, while an exhaustive 
description of the tool will be given in a forthcoming paper 
(Fasano et al., in preparation). Here we just mention that, among 
the 14 newly devised diagnostics, the most effective one in order 
to disentangle ellipticals from SO galaxies turned out to be an 
Azimuthal coefficient, measuring the correlation between azimuth 
and pixel flux relative to the average flux value of the elliptical 
isophote passing through the pixel itself. Figure [3] illustrates 
the distributions of the Azimuthal coefficient for the samples of 
ellipticals and SO galaxies used to calibrate MORPHOT. 

More importantly for our purposes, the quantitative discrep- 
ancy between automatic (MORPHOT) and visual classifications 
turns out to be similar to the typical discrepancy among visual 
classifications given by experienced, independent human classifiers 
(r.m.s.~1.5-2.5 T types). This was extensively tested on a calibra- 
tion sample of 931 WINGS galaxies with both visual and MOR- 
PHOT classification, where the visual classification was carried out 
independently by GF and AD. This is illustrated in Figure [4] where 
the automatic classification is also shown to be bias-free in the over- 
all range of morphological types. 

For now, we can apply MORPHOT just to the WINGS imag- 
ing, because the tool is calibrated on the WINGS imaging charac- 
teristics, and we defer to a later time a more generally usable ver- 
sion of the tool. To verify directly that the two methods adopted at 
different redshifts (see Section l4~2l are consistent, we can apply the 
same "method" (visual classification and persons) that was used at 
high-z on the low-z images. 

To this aim, 3 of the classifiers that in 2007 visually classified 
all the EDisCS galaxies (BMP, AAS, VD) now performed a visual 
classification of WINGS galaxies. This was done on the subset of 
WINGS galaxies that was used to calibrate MORPHOT on the vi- 
sual WINGS morphologies, including only galaxies that enter the 
sample we analyze in this paper (173 galaxies). 

The results (see FigO show agreement between the three 



broad morphological classes assigned by the EDisCS classifiers 
with the WINGS visual classification in ~83% of the cases, and 
with MORPHOT in ~75% of the cases. 

Again, these discrepancies turn out to be similar to the typi- 
cal discrepancy among visual classifications given by experienced, 
independent human classifiers, so we conclude that the different 
methods adopted provide a comparable classification. 



4.1.1 Mass limited sample 

The spectroscopic magnitude limit of the WINGS survey is V=20. 
Considering the distance module of ~ 37.5 of the most distant 
WINGS cluster, and a reddest color of (B — V) = 1.2, the magni- 
tude limit corresponds to a mass limit M» = 10 9 8 Mq, above which 
the sample is unbiased. Adopting this limit, the final sample com- 
posed by galaxy cluster members, within 0.6i?200, excluded the 
BCGs, consists of 1233 galaxies, of which 396 ellipticals galax- 
ies, 586 SOs, 239 late type galaxies and 12 galaxies with no clas- 
sification (see Table 3). The corresponding numbers weighted for 
completeness are also given in Table 3, for a total of 1894.42 galax- 
ies. Numbers of WINGS galaxies above the EDisCS mass limit 
(M, = 10 10 ' 2 M©, see below) are also listed in Table 3. 



4.1.2 Magnitude limited sample 

Since we wil l use the morphological fractions derived by 
IPoggianti et all (2009), we use the magnitude limit adopted in that 
study, which is M(/=-19.5, corresponding to the passively evolv- 
ing limit of the E DisCS morphological study at higher redshifts 
jPesai et al]|2007|) . The final sample of members within 0.6^200, 
excluded the BCGs, consists of 263 ellipticals, 410 SOs, 191 late- 
type galaxies, and 10 galaxies with no morphological classification, 
for a total of 874 galaxies, as listed in Table 3. Numbers weighted 
for completeness, for a total of 1340.73 galaxies, can be found in 
Table 3. 
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Figure 5. Comparison between the visual classification performed by the EDisCS classifiers and the WINGS visual classification (left panel) and the automatic 
classification performed by MORPHOT (right panel). The circle radius is proportional to the number of galaxies. 



4.2 High-z sample: EDisCS 

For EDisCS, there are 605 spectroscopically confirmed members. 
381 of them (~ 63%) have HST morphologies. The spectroscopic 
magnitude limit ranges between 1=22 and 1=23 depending on red- 
shift, and the most conservative spectroscopi c mass limit using a 
iKroupal feOOll) IMF becomes M = 10 10 6 M Q ^Vulcani et al1l20K]h . 
Therefore, relying only on spectroscopy, the lower mass limit of our 
study would be very high, and numbers would be relatively small. 

For this reason, we decided to use all photo-z members (see 
92.2| l. All the EDisCS mass functions analyzed below are based 
solely on the photo-z membership, unless otherwise stated. The 
photo-z technique allows us to push the mass limit to much lower 
values than the spectroscopy (see below). As shown in Appendix A, 
spectroscopic and photo-z techniques give very consistent results in 
the mass range in common. 

As for the WINGS sample, BCGs and galaxies at radii greater 
than r = 0.6^200 have been excluded from the analysis. 

Morphologies are discussed in detail in iDesai et al. I fe007t) . 
The morphological classification of galaxies is based on the visual 
classification of HST/ACS F814W images sampling the rest-frame 
~ 4300-5500 A range, similarly to WINGS. 

4.2.1 Mass limited sample 

The magnitude completeness limit of the EDisCS photometry is 
/ ~ 24, though the completeness cle arly remains high to magni- 
tudes significantly fainter than / = 24 JWhite et alj|2005h . 

We consider the most distant cluster, cl 1216.8-1201 , that is 
located at z ~ 0.8 and determine the value of the mass of a galaxy 
with an absolute B magnitude corresponding to / = 24, and a color 
(B — V) ~ 0.9, which is the reddest color of galaxies in this cluster. 
In this way, the EDisCS mass completeness limit based on photo-z 
isM» = 1O 1O2 M©0 

The completeness is confirmed by the analysis of 

3 Note that with this selection only a few (1.7%) galaxies in the mass- 



Rudn ick et all (2009), who find that for the magnitudes used 
in this paper (M g always ^ -20.2) the photo-z counts are fully 
consistent with the statistically background subtracted counts for 
both red and blue galaxies. 

The final mass-limited EDisCS sample above this limit con- 
sists of 489 galaxies, 156 of which are classified as ellipticals, 64 
as SOs, 252 as late-type galaxies, and 17 galaxies with unknown 
morphology (see Tabled- 

4.2.2 Magnitude limited sample 

We us e the same magn itude limit as in the EDisCS morphological 
study (Des ai et al . 2007), My =-20, which for passive evolution cor- 
responds to the magnit ude limit in the WINGS morphological study 
( tPoggianti et al.| [2009 ) . With this limit, our sample consists of 544 
galaxies, 151 of which are ellipticals, 67 SOs, 304 late-types, and 
22 galaxies with unknown morphologies (see Table[3). 



5 RESULTS: THE MASS FUNCTION 

In this section we analyze only the mass limited samples of both 
data-setsQ 

We use galaxy stellar masses to characterize the mass distri- 
bution of different galaxy types in clusters, both at low- and high-z. 

Above the completeness limit, we build histograms to define 
the mass distribution. In each mass bin, we sum all galaxies of all 
clusters to obtain the total number of galaxies, then we divide this 
number for the width of the bin, to have the number of galaxies per 
unit of mass. The width of each mass bin is 0.1 dex. In the case 
of spectroscopic data, as for WINGS, in building histograms each 

limited sample are fainter than 1=23, and therefore have no morphological 
classification from Desai et al. 2007. 

4 Appendix B presents the mass functions for magnitude limited samples 
of WINGS and EDisCS galaxies, and clearly shows the biases inherent in 
this type of selection. 
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Table 3. Number of galaxies above the different limits. For WINGS both observed numbers and numbers weighted for spectroscopic incompleteness are given 



galaxy is weighted by its incompl eteness correct ion. Errorbars are 
computed using poissonian errors (Gehrels 1986). 

Clusters of different masses could have a different role in the 
total mass function, so we have tried to test the importance of this 
effect. For WINGS galaxies, we have determined the mass func- 
tions separately for clusters with a < 700 kms~\ 700 kms~ l < a < 
800 kms~ and a > 800 kms~ l (plots not shown) and we found no 
differences. Similarly, for EDisCS galaxies we have tried building 
the total mass function normalizing for the number of galaxies in 
each cluster, to avoid the possible larger contribution of larger clus- 
ters to the total mass function. Once again, we have found (plots 
not shown) no differences with the mass function presented in this 
paper. So, a different combination of the mass functions of different 
clusters could only influence the final normalization; anyway, since 
we are interested mainly in the shape of the mass distributions and 
not in the counts (see below), this doesn't affect our results. Con- 
sequently, we conclude that the mass of the cluster, over the cluster 
mass ranges considered here, does not influence significantly the 
shape of the total mass function, and in the following we use the 
simple sum of all galaxies in all clusters. 

To quantify the differences between different mass functions, 
we perform Kolmogorov-Smirnov (K-S) tests. However, the stan- 
dard K-S test does not consider the completeness issues when it 
builds the cumulative distribution (since it assigns to each object 
a weight equal to 1). So we modified the test, to make the relative 
importance of each galaxy in the cumulative distribution depend on 
its weight. In the following, we will always use this modified K-S 
test. Obviously, for photo-z data all galaxies have a weight equal 
to 1 , and using the modified test is equivalent to using the standard 
one. 



5.1 The total mass function and its evolution 

We first compare the mass distribution of all galaxies at low and 
high redshift. As we are interested in the shapes of the mass func- 
tions and not in the total number of galaxies, we normalize the 
mass distribution of EDisCS galaxies to the mass distribution of 
WINGS galaxies to have the same total number of galaxies above 
log M„/Mq ^ 11. In this way we assume th at the most massive 
galaxies are already in place at z~ 1 ( see, e.g jFontana et alj2004l ; 
IPozzetti et aT]|2007l) . 

In Figure|6]the mass distribution of WINGS galaxies (crosses) 
is compared to the mass distribution of EDisCS galaxies (squares). 
The comparison is meaningful only above the highest of the two 
mass completeness limits, that is for logM*/M ^ 10.2. 

The two mass distributions are very different. Indeed, the 
EDisCS mass function is flat for logM»/M s: 10.8. The WINGS 
mass function is rather flat up to about logM»/M ~ 10.6 and then 
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Figure 6. Comparison of the mass distribution of EDisCS (blue filled 
squares) and WINGS (black crosses) galaxies of all morphological types, 
for the mass limited sample. Mass distributions are normalized to have the 
same total number of objects at log M»/M ^11. Errors are poissonian. 
Only points above the respective completeness limits are drawn. Numbers 
in the labels are the total number of galaxies above the highest of the two 
completeness limits. For WINGS, both observed (first) and weighted (sec- 
ond) numbers are given. The inset shows the cumulative distributions, with 
the relative K-S probability, for galaxies above the common mass limit. 
The mass distributions depend on redshift: for a similar number of higher 
mass galaxies, for logM„/M ^ 10.8 there are proportionally less massive 
galaxies at low-z than at high-z. 



it begins to decline more steeply, with a massive end slope similar 
to that of the EDisCS mass function. The striking difference is the 
dearth of galaxies less massive than logM»/M ^ 10.8 in EDisCS 
compared to WINGS, for the same number of massive galaxies. 
The K-S test confirms that the two distributions are different at the 
99.99% level. We remind that the results of the K-S test do not de- 
pend on the normalization we adopted in the plots. So, also using 
other normalization criteria, we would reach the same conclusions. 

We conclude that the total mass distribution of galaxies in 
clusters evolves with redshift. In a mass-limited sample, cluster 
galaxies at high-z are on average more massive than in local clus- 
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ters. Assuming no evolution at the massive end, the population of 
logM, /Mq ^10.8 galaxies must have grown significantly between 
z = 0.8 andz = 0. 

In principle, the evolution of the galaxy stellar mass function 
can be due to several factors: a) strong environmental mass segre- 
gation, for which galaxies infalling into clusters between the two 
redshifts have a different mass distribution than galaxies already 
in clusters at high-z (infalling galaxies being on average less mas- 
sive); b) mass loss from massive galaxies due to harassment; c) 
galaxy merging; d) mass growth due to star formation. In §7 we 
will discuss which effects are likely to play an important role and 
which ones can probably be neglected in this study. 



5.2 The mass functions of different galaxy types 

We now analyze separately the mass function of each morpholog- 
ical type, to see if and how the distribution of masses depends on 
the galaxy type. In this case, we do not apply any normalization 
among the different mass functions, to show what morphological 
type dominates as a function of mass. Figure[7]shows the mass dis- 
tributions for galaxies of the WINGS dataset (left panel) and of the 
EDisCS dataset (right panel). Galaxies of different morphological 
types are indicated by different colors and symbols. 

In Table [4] we give the analytic expression for the mass func- 
tions of the differ ent morphologica l types for both samples, esti- 
mating the best-fit lSchechtea d 19761) parameters (a, M* , </)*)■ 

For WINGS, the shape of the mass distribution strongly de- 
pends on morphological type. The total mass function is domi- 
nated by SO galaxies up to logM./M© ~ 11, and by elliptical 
galaxies at higher masses. The total mass function is quite flat up 
to logM„/M0 ~ 10.6, as it is the function of both SOs and el- 
lipticals. Instead, the mass function of late-type galaxies declines 
over the whole mass range. While there are elliptical galaxies 
at all masses, there are no SOs or late-types more massive than 
logM»/MQ ~ 11.4 and logM,/M Q ~ 1 1 .6 respectively. 

We perform a K-S test, both for a mass limit logM, /Mq 9.8 
(the WINGS limit), and for logM*/M Q > 10.2 (the same limit as 
EDisCS). The K-S results are summarized in Table [5p| 

For logM„/M0 ^ 9.8, the morphological mass functions sig- 
nificantly differ one from another: ellipticals, SO and late-type 
galaxies all have different mass functions, with the K-S that ex- 
cludes the possibility that they are drawn from the same par- 
ent distribution with a probability always greater than 99.4%. For 
logM»/Mo 10.2, SOs and late-type galaxies are not statistically 
distinguishable any more (probability = 91.4%). This suggests that, 
except for the least massive galaxies, the SO and spiral mass func- 
tions could in fact be similar. This could be due to the fact that 
when we analyze the sample at logM^/M© 10.2 we are actually 
excluding Sc and Irregulars galaxies, that mostly have masses be- 
low this limit. However, in our sample, they represent only 5% of 
the total population of late-type galaxies, so we do not expect they 
matter a lot. 

Comparing each type with all galaxies, for logM*/Mo ^ 9.8 
only SO galaxies could have a similar distribution to the sum of all 



5 We note that with the term "all" we refer to all galaxies regardless of their 
morphological type and we are including also galaxies for which we do not 
have a morphological classification, that are only a few (see Table 3). How- 
ever, we have checked that the mass distribution is the same considering all 
galaxies or only the sum of ellipticals, SOs and late-type galaxies. 
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Figure 8. Mass distribution of EDisCS late galaxies split into Sa+Sb (black 
crosses) and Sc+lrr (blue filled squares), for the mass limited sample. Mass 
distributions are normalized to have the same total number of objects at 
logM, /Mq ^11. Errors are poissonian. Numbers in the labels are the total 
number of galaxies above the completeness limits. The inset shows the cu- 
mulative distributions, with the relative K-S probability, above the common 
mass completeness limit. The two distributions are quite different, indicat- 
ing that early- and late-spirals have different mass distributions. 



galaxies (Pks ~ 70%), while in the other two cases the test indi- 
cates that the distributions are different, with a probability at least 
> 98.5%. Instead, above logM„/M Q = 10.2 the K-S test is less 
conclusive: it rejects the hypothesis of similarity only between el- 
lipticals with all galaxies, with a probability of 99.9%. This demon- 
strates how important it is to reach deep mass limits: in fact the 
greatest differences among the mass functions tend to emerge at 
lower masses. 

For EDisCS (right panel in Figure the total mass func- 
tion is dominated by the mass function of late-type galaxies up 
to logM»/M < 11 (and not by SOs galaxies as in WINGS). The 
most massive bins of the total function are dominated by elliptical 
galaxies. In fact, there are no SOs or late-types more massive than 
logMj/M© ~ 11.5 andlogM./M© ~ 1 1.6 respectively. The mass 
function of SOs is very noisy, due to their low number in the sample. 

At high-z, the data cannot rule out that the mass functions of 
the different morphological types are all similar, with K-S proba- 
bilities always smaller than or at most equal to 90% (see Tabled- 
We also split late-type galaxies in Sa+Sb and Sc+lrr galaxies, to in- 
vestigate if the two populations behave differently (Figure[8}. Even 
though the K-S seems to be not conclusive (it gives a probability of 
~ 93%), we find that mass functions are quite different (Figure[8]l, 
indicating that there are proportionally more lower mass galaxies 
(logMt/M© < 1 1) among Sc+lrr than among Sa+Sb| 

Also the comparison between each morphological type and all 



We note that such an analysis cannot be performed in WINGS, due to the 
very small number of Sc+lrr. 



WINGS 



Galaxy stellar mass functions and morphologies 1 1 

EDisCS 



-I 1 r- 



T" 



T" 




xoll (1233, 1894.418) 
• ellipticals (396, 603.314; 

•irSO (586, 912.452) 
Olate (239, 358.516) 



t — I — i — i — i — r 




1 1.5 



o 



' • o x o x 




xoll (489) 

• ellipticals (156) 

irSO (64) 
Olate (252) 



ON 



-- If-, if : - 
U 



;>.. 6-,- 



all-ell 62.9% 
all-SO 82.9% 
all-late 14.8% 
ell-SO 53.2% 
ell-late 90.3% 
SO-late 80.4% 

_l i i 



<) T 



10.5 



1 1 

Log, M 



1 1.5 



Figure 7. Mass distribution of WINGS galaxies (left panel) and EDisCS galaxies (right panel). In the insets, the cumulative distributions are shown. For 
WINGS, the right inset is done for the EDisCS mass limit, logM, /Mq = 10.2. Black crosses represent all galaxies, filled red circles elliptical galaxies, 
green stars SOs and blue open circles late-type galaxies. Errors are defined as poissonian errors. Numbers in the labels are the number of galaxies in each 
morphological class, above the res pective mass limit . For WINGS, both observed (first) and weighted (second) numbers are given. For EDisCS, the K-S 
probabilities are also given. Best-fit Schechter 1 1976) mass function are showed dotted. The values of the parameters of the fit are presented in Table [4] At 
low-z, different morphological types have different mass functions, while at high-z no conclusion can be drawn. 







logM* 


a 


0* 


WINGS 


all galaxies 


11.667±0.052 


-0.987±0.009 


1.552±0.055 




early-types 


11.643±0.055 


-0.977±0.011 


1.526±0.059 




ellipticals 


11.879±0.076 


-1.017±0.014 


1.13±0.067 




SOs 


11.391±0.052 


-0.925±0.020 


1.680±0.089 




late-types 


11.683±0.171 


-1.04±0.031 


0.992±0.138 


EDisCS 


all galaxies 
early-types 
ellipticals 
SOs 

late-types 


11.684±0.061 
11.700±0.081 
11.847±0.126 
11.300±0.125 
11.523±0.067 


-0.915±0.026 
-0.896±0.041 
-0.947±0.044 
-0.662±0.147 
-0.875±0.044 


1.577±0.105 
1.423±0.138 
1.159±0.148 
1.679±0.235 
1.560±0.141 



Table 4. Best-fit Schechter ( 1976) parameters (a, M* , </>*) for the galaxy mass functions of the different morphological types discussed in this paper. Fits are 
perfomed directly on the logarithmic points. 



WINGS 



EDisCS 





P K - S {M* 10"- 8 M Q ) 


P K - S (M* 10 la2 M Q ) 


P K -s(M, > 1O 1O2 M ) 


ell-SO 


99.84% 


99.99% 


53.18% 


ell-late 


99.44% 


98.20% 


90.20% 


SO-late 


99.72% 


91.40% 


80.43% 


all-ell 


99.38% 


99.9% 


62.98%% 


all-SO 


69.52% 


95.53% 


82.89% 


all-late 


98.64% 


24.94% 


14.77% 



Table 5. Results of the Kolmogorov-Smirnov test performed on the mass function of different morphological types for the WINGS and EDisCS mass limited 
sample. For WINGS, results also above the EDisCS mass completeness limit are given. Pk-s is the probability that two distributions are drawn from different 
parent distributions. The K-S has been modified to take into account the different weights of galaxies observed spectroscopically. 
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galaxies does not give conclusive results, because the K-S probabil- 
ities are always smaller than 90%. In WINGS, we have seen that, 
for the same mass limit, only SOs and late-types could have the 
same mass distribution, while they both differ from ellipticals. The 
lack of significance of the differences in EDisCS might be due to 
poorer number statistics, as well as to the higher mass limit. 

To check if the lack of significance in EDisCS is real or 
above all due to poor statistics, we performed 100 Monte Carlo 
simulations extracting randomly from the WINGS sample at 
logM*/MQ ^ 10.2 a subsample with the same number of galax- 
ies as EDisCS (that is 156 ellipticals, 64 SOs and 252 spirals). We 
found that the probability of the K-S test to be conclusive is very 
low in all cases (at most in 56% of the simulations when we com- 
pare ellipticals and late-types). These results are consistent with 
the hypothesis that in EDisCS the lack of a significant difference 
between the mass functions of different types could be due to poor 
number statistics. 

The conclusion we can draw from this section is that, at low 
redshift, different morphological types have significantly different 
mass functions, although SOs and late-type galaxies more massive 
than logM*/M0 = 10.2 may have similar mass distributions. No 
conclusion can be drawn at high redshift. 

5.3 Evolution of the mass functions of different 
morphological types 

We now directly compare the mass distribution of each morpholog- 
ical type in the two reshift bins, with the aim to understand if the 
shape of the mass function of each given type evolves with redshift, 
and what type of galaxies drives the evolution of the total mass 
function. As for the total mass function, we normalize the mass 
distribution of EDisCS galaxies to the mass distribution of WINGS 
galaxies to have the same total number above log M*/Mq ^11, 
since we want to compare the shape of the distributions. 

Figure [9] compares the mass distribution of galaxies of each 
morphological type at high and low redshift. The mass distribution 
of ellipticals is shown in the top left panel. For both WINGS and 
EDisCS, the mass distribution starts rather flat at low masses, with 
a hint for a possible dip in EDisCS, then begins to steepen at higher 
masses. Statistically, the EDisCS and WINGS distributions are sig- 
nificantly different. As for the total mass function, there is a deficit 
of less massive galaxies (logM, /Mq ^ 10.6) at high-z compared 
to low-z, for a given number of more massive galaxies. A K-S test 
rejects the null hypothesis of similarity of the two distributions with 
a probability >99.9%. 

In the top right panel of Figure [9] the mass distribution of SO 
galaxies is shown. It changes dramatically with redshift. WINGS 
galaxies have a flat mass function up to log M*/Mq ~ 10.6 and 
then this becomes steeper. Instead, the mass function of EDisCS 
galaxies shows a rise, a peak at about logM, /Mq ~ 10.6, and a fall. 
The K-S test states that the mass distributions at the two redshifts 
are driven from two different parent distributions with a probability 
> 99.99%. This result emerges with a high statistical significance 
even though SO galaxies are very rare in the EDisCS sample. 

In the bottom left panel of Figure|9]we analyze the mass func- 
tion of early-type galaxies. This can be directly compared with the 
mass function of late-types (bottom right panel). The differences 
existing in the total mass distribution still remain evident in the 
mass function of early-type galaxies. While the mass distribution of 
high mass galaxies is rather similar at the two redshifts, the dearth 
of galaxies with logM, /Mq ^ 10.9 in EDisCS is even more evi- 
dent than in the total mass function. Once again, the K-S test rejects 



the possibility that the two mass functions derive from the same dis- 
tribution with a probability >99.9%. Of course, this behavior was 
expected, being the early-type galaxies simply the sum of ellitpicals 
and SOs. 

The bottom right panel of Figure[9]shows the mass distribution 
of late-type galaxies. The shape of the mass distribution of EDisCS 
and WINGS is different up to logM, /Mq ~ 1 1 and then it is rather 
similar at higher masses. For WINGS late-types, the mass distribu- 
tion is immediately declining, while that of EDisCS starts flat. The 
K-S test gives a probability of ~ 98% that the two distributions are 
not driven from the same parent distribution. 

The conclusion we can draw from this section is that the mass 
distribution of each morphological type evolves with redshift. All 
types have proportionally more massive galaxies at high- than at 
low-z. 

In addition to the processes mentioned in §5.1 (infall, harass- 
ment, merging and star formation) that can change the total mass 
function, there is an additional effect that could influence the ob- 
served evolution of the mass functions of each morphological type: 
the morphological transformation from one type to the other. 



6 RESULTS: THE MORPHOLOGICAL FRACTIONS 
6.1 The Morphology-Mass relation 

It is interesting to analyze how the incidence of each morphological 
type (the number of galaxies of that type over the total) changes 
as a function of the stellar mass (hereafter the morphology-mass 
relation, MMR). 

As in the following w e will c omp are our findings with the re- 
sults found bv lDesai etall d2007l) and IPoggianti et alj J2009I) that 
use magnitude limited samples, we discuss the MMR both in mass 
limited and in magnitude limited samples. Since above the com- 
pleteness limit of our samples the fraction of objects without mor- 
phological classification is negligible (at most 4%) in this analysis 
we do not consider them, so we refer to "all" galaxies simply as the 
sum of early and late- type galaxies. 

Figure [10] shows the fractions of the different morphological 
types as a function of galaxy mass, for both EDisCS and WINGS. 
In the left panel, the mass li mited samples are presented. Errors are 
defined as binomial errors lGehrelslll986l) . 

The trend of ellipticals at high and low redshift is similar: the 
elliptical fraction is flat at low masses, then it increases at higher 
masses. Ellipticals represent more than 50% of galaxies at any mass 
greater than logM„/M ~ 11 at low-z and logM./M© ~ 11.4 at 
high-z. Thus, WINGS elliptical galaxies start to dominate at a lower 
mass than in EDisCS. The most massive galaxies tend to be almost 
all ellipticals, at both redshifts. 

The morphology-mass relation of SOs and late-types galax- 
ies depends even more strongly on redshift. In the Local Universe, 
for masses up to logM„/M0 ~ 11 SOs are by far the most com- 
mon galaxies, and their fraction is rather constant with mass up to 
logM„/M0 ~ 10.9, whereas it decreases steeply towards higher 
masses. In contrast with the local Universe where SOs are much 
more frequent at low than at high masses, at high-z SOs are equally 
rare (^ 20%) at all masses. 

Late-type galaxies are the least represented morphological 
class in the Local Universe (^ 20% at all masses), their trend 
hardly depending on mass (excluding a possible fluctuation at 
logM„/M Q ~ 10.9). At higher redshift, the situation is differ- 
ent: late-types are the most common type of galaxies up to 
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Figure 9. Comparison of the mass distribution of EDisCS (blue filled squares) and WINGS (black crosses) galaxies, for the different morphological types. 
Mass distributions are normalized to the WINGS total number of objects with logM„/Mo ^11. Errors are defined as poissonian errors. Numbers in the labels 
are the number of galaxies in each morphological class, above the respective mass limit. For WINGS, both observed (first) and weighted (second) numbers 
are given. The cumulative distributions are shown in the insets. Each morphological type has a mass distribution that depends on redshift. In WINGS there are 
always proportionally more less-massive galaxies than in EDisCS. 



logM, /Mq ~ 11.3, with no clear trend with mass at lower masses, 
and a steep decline at higher masses. 

In magnitude limited samples, the dependence of the morpho- 
logical fractions on mass is even more pronounced, at both red- 
shifts. At low masses, late-type galaxies dominate the total pop- 
ulation, while the presence of elliptical galaxies is negligible. On 
the contrary, at higher masses elliptical galaxies are the most nu- 



merous, and late-type galaxies the least numerous type of galax- 
ies. As far as SOs are concerned, while in WINGS they domi- 
nate at intermediate masses (logM»/M Q ~ 10.1-10.9), in EDisCS 
they do not dominate at any mass. Comparing high and low red- 
shift, in general there are many more intermediate-mass and mas- 
sive EDisCS late-types than WINGS late-types, that instead dom- 
inate only the first bins up to logM„/M0 ~ 10.1, and then fall to 
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Figure 10. Morphological fractions for the mass limited (left panel) and magnitude limited (right panel) samples for galaxies at different redshifts. Solid 
symbols and lines refer to EDisCS, empty symbols and dotted lines refer to WINGS. Squares refer to ellipticals, circles refer to SOs, triangles refer to late- 
types. Errors are defined as binomial errors. In both samples, the morphological mix changes with galaxy stellar mass and with redshift at a given mass (see 
text for details). 



zero by logM./Mo ~ 11.5. Similarly to the mass-selected sam- 
ple, WINGS ellipticals start to dominate at lower masses than in 
EDisCS. 

These results show how the proportion of galaxies of differ- 
ent morphological types changes with galaxy mass. This MMR is 
significantly different in our mass-limited and magnitude-limited 
samples. In the mass-limited sample, the incidence of each type 
does not depend on mass between our completeness mass limits 
and logAf»/AfQ ~ 11-11.1 (flat trend), while it strongly changes 
at higher masses. Both for the mass-limited and the magnitude- 
limited samples, the MMR evolves strongly with redshift, in con- 
nection to the evolution of the morphological fractions described 
below. 



6.2 Evolution of morphological fractions 

The morphological fractions for both mass and magnitude limited 
samples are given in Table|6] WINGS fractions are given also above 
the EDisCS mass limit. For WINGS, fractions are computed both 
for unweighted and for completeness-weighted numbers. We can 
note that the completeness correction does not change the resulting 
fractions significantly. 

Except for late-types, in WINGS, fractions of different mor- 
phological types do not depend on the choice of the mass limit, 
remaining constant within the errors. Moreover, the choice of the 
sample (mass or magnitude limited) does not alter the trends ob- 
served. 

The fraction of elliptical galaxies does not change with red- 
shift, being almost constant (~ 30%). Instead, the fraction of SO 
and late-type galaxies considerably changes, being SOs much more 
common in the Local Universe, while spirals are proportionally less 
common. 





WINGS 


EDisCS 




My C -19.5 


M v -20 




(Poeeianti et al. 2009) 


(Desai et al. 2007) 


ellipticals 


33 ± 7% 


29 ± 1% 


SOs 


44 ± 10% 


14 ± 1% 


early-types 


77 ± 12% 


43 ± 2% 


late-types 


23 ± 9% 


57 ± 3% 



Table 7. Morpho logical fractions obtained from Pogg ianti et aD [2009); 
lDesaietalJj2007h . 



As a consequence, the early-type fractions strongly evolves, 
almost doubling from high- to l ow-z. This is in agreement with 
the finding o f Kovac et al for the field, but it is in contrast 

with lHolden et alj d2007r T who find no evolution in the early-type 
fraction in clusters from z ~ 0.8 to the current epoch in a mass 
limited sample with logM./M© 10.55. Our trend still remains 
outstandin g ( ~ 50% at z~ .8 and ~ 83% at z~ ) even if we 
adopt the lHolden et al.1 d2007l) mass limit. 

Our r esults are in very good agree ment with those fo und pre- 
viously by IPoggianti etafl d2009l) and iDesai et al.l d2007t) for the 
WINGS and EDisCS datasets, respectively, listed in Table|7] 

We note that here we use only spectroscopic member s in a 
subset of 21 WINGS clusters, while IPoggianti et alj ( f2009h used 
all galaxies in the photometric catalog of 72 WINGS clusters. 
Moreover, we note that considering EDisCS photo-z m embers, our 
fractio n are in perfect agreement with those found by Des ai et al.l 
d2007l) . who used the total photometric catalog and a statistical 
background-subtraction technique, indicating that photo-z tech- 
niques to assign membership and statistical subtraction give con- 
sistent results. 
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M„ > 10 




WINGS 
M, 5: 1O 1O - 2 M 




-19.5 


EDisCS 

M» > 1O 1O ' 2 M 


My ^ -20 




%obs 








%obs 




% 


% 


ellipticals 
SOs 
early 

late-types 


32.4±1.3% 
48.0±1.4% 
80.4±1.1% 
19.6±1.1% 


32.2±1.1% 
48.7±1.2% 
80.9±0.9% 
19.1±0.9% 


33.2±1.7% 
50.5±1.9% 
83.7± 1.4% 
16.3±1.4% 


32.7±1.4% 
51.5±1.5% 
84.2±1.1% 
15.8±1.1% 


30.4±1.6% 
47.4±1.7% 
77.8±1.4% 
22.2±1.4% 


30.0±1.3% 
48.4±1.4% 
78.4±1.1% 
21.6±1.1% 


33.0±2.2% 
13.6±1.6% 
46.6±2.3% 
53.4±2.3% 


28.9±2.0% 
12.8±1.5% 
41.8±2.1% 
58.2±2.1% 



Table 6. Morphological fractions of galaxies in both mass and mag-limited samples. Errors are computed as binomial errors. For WINGS, both observed and 
completeness-weighted numbers are listed. Since galaxies with unknown morphology are always ^ 4%, they are not included in the computation of fractions. 





WINGS 


EDisCS 




M, 10' 8 M Q 


M, ^ 1O 1O ' 2 M 


M» 3= 10 10 ' 2 M Q 


ellipticals 


42 ± 4% 


43± 5% 


39 ± 5% 


SOs 


42 ± 4% 


42 ± 4% 


13 ± 2% 


early-types 


84 ± 6% 


85 ± 7% 


52 ± 5% 


late-types 


16 ± 2% 


15 ± 2% 


46 ± 4% 


unknown 


0.63 ± 0.22% 


0.54 ± 0.24% 


2± 1 % 



Table 8. Fraction of mass in galaxies of each type, above the different lim- 
its for WINGS and EDisCS. Errors are calculated with the propagation of 
errors, assuming a mean error on each galaxy mass equal to the mass itself. 



6.3 In what type of galaxies is most of the mass? 

In Table [8] we present the amount of mass in each morphological 
type, to quantify the relative contribution of each class to the total 
stellar galaxy mass. We stress that only galaxies more massive than 
our mass limits are considered, and no attempt is made to extrapo- 
late to lower galaxy masses to obtain the total mass, integrated over 
all galaxy masses. 

Errors are calculated with the propagation of errors, assuming 
a mean error on each galaxy mass equal to the mass itself, which 
corresponds to the 0.3 dex error in mass discussed in §3. 

At low redshift, for logM, /M Q ^ 10.2, ~ 85% of the mass is 
in early-type galaxies, of which half in SOs and half in ellipticals. At 
high redshift, only about 50% of the mass is in early-type galaxies, 
of which most (~ 3/4) is in ellipticals. 

Interestingly, the fraction of stellar mass in WINGS SOs is the 
same as in EDisCS late-types. So it is not just the number frac- 
tion of high-z spirals (SOs) that matches that of low-z SOs (spirals), 
but also the mass fractions. Consequently, it is not just the number 
fraction of ellipticals that remains constant with redshift (see previ- 
ous section), but also the fraction of mass in ellipticals, being about 
40% both in EDisCS and in WINGS. 



7 DISCUSSION 

In this paper, we have presented mass functions both of all galaxies 
and of different morphological types in clusters at different epochs. 

For a high mass threshold (logM„/MQ ~ 11), the shapes of 
our z = 0.7 and z = 0.05 mass functions are similar, consistent with 
high-mass galaxies having already established their mass at high- 
z. The K-S test supports this assumption, giving a probability ~ 
90% that the two distributions are not driven from the same parent 
distribution. 

This is in ag r eemen t wit h results for the field, whe re for exam- 
ple, iBundv et all J2006h and iFranceschini et alj fcOOfJ) found that 



the high-mass end of the mass function is quite constant at least 
from z ~ 1 to z ~ 0. This is also in agreement with previous studies 
of the bright end of the galaxy luminosity function in clusters at 
red wavelengths (De Propris et al. 2007, Toft et al. 2004, Andreon 
2006). 

In contrast, going to less massive galaxies, the differences 
increase: clusters in the local Universe have proportionally many 
more less-massive galaxies than clusters in the distant Universe. 

We have found the same when the sample is divided by mor- 
phological type. Furthermore, we have shown how the morpholog- 
ical mix changes with galaxy stellar mass, and with redshift at a 
given mass. 

Both the number fraction and the mass fraction of ellipticals 
does not change with z, while the fractions of SOs and spirals are 
inverted at the two redshifts. 

At the end of §5.1 and 5.3 we have identified some mecha- 
nisms that in principle could contribute to the observed evolution 
of the mass functions and of the morphological mix. 

In this section, we discuss each one, trying to understand 
which ones can be the most relevant in driving the trends observed. 

• merging: a galaxy can gain mass through mergers with other 
galaxies; so mergers could drive both galaxy mass assembly and 
dynamical evolution in a hierarchical scenario. However, in clus- 
ters this process is not expected to play an impor tant role (see e.g. 
lOstrike jfl980l : iMakino & Hutll 19971 ; lMihosll2004l) : th e high veloc- 
ity di spersions in clusters suppress mergers (see e.g. lEllison et al.l 
2010 and references therein). More importantly, this type of process 
should favor the growth of the high-mass end of the mass function 
relative to the low-mass end, as we observe. Finally, visual morpho- 
logical classificatio ns show that the fr action of colliding galaxies in 
clusters is low (e.g. lDesai et alj2007l) . 

• harassment: galaxies in clusters are subject to frequent high 
speed encounters that strip a galaxy of part of its mass and drive 
a morphological transformation. Harassment has the potential to 
change any internal property of a galaxy within a cluster including 
the gas distributio n and content, the o rbital distribution of stars and 
the overall shape l lMoore et alJI 1996b . Its strength depends on the 
collisional frequency, on the strength of the individual collisions, 
on the cluster's tidal field and on the distribution of the potential 
within galaxies. It has long timescales since it needs multiple en- 
counters to be efficient in removing gas and stars and quenching 
star formation. It effectively perturbs low-mass gala xies while its 
effects on massive objects sho uld be less pronounced ( Moo re et"al] 
ll996l ; lBoselli & Gavazzill 20061) . 

• environmental mass segregation of in] "ailing galaxies: clusters 
are far from being closed boxes, so the infall of galaxies from sur- 
rounding areas is a mechanism that certainly takes place. If galax- 
ies surrounding the clusters have a steeper mass function than the 
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cluster, they would make a greater contribution to the total mass 
function at intermediate to low-masses. 

Characterizing this process is very complex, since there are sev- 
eral factors that must be understood, a) First of all, we should 
know the rate of infalling galaxies with time, to understand how 
the number of galaxies in clusters increases with redshift. b) Then, 
we should know the relative fraction of the different morpho- 
logical types of infalling galaxies. It depends on the morpholog- 
ical mix of galaxies located in the surrounding areas the clus- 
ters. In the field th e morphological mix evolves with redshift 
dPannella et al . 2006), so we also need the proportions of the mor- 
phological t ypes of infalling galaxie s to change from z ~0.8 to z 
~0 (though Ivan der Wei et alT ( l2007l) find that the early-type frac- 
tion in the field and group environments is quite constant (~45%) 
from z ~ 0.8 to the current epoch). In the outskirts of the cluster, 
the gl obal morphological mixing i s correlated with local density 
(e.g. iHuertas-Companv et alj|2009l) . In general, in the field there 
are more late-type galaxies than early-types (at z ~ 0.6, ~ 70% at 
low masses, ^ 50% at intermediate masses (Oesch et al. 2 oToh ) so 
late type galaxies are probably the favourite infalling galaxies, c) 
Most importantly, we should know the exact mass distribution of 
infalling galaxies, that c an be expected to dif fer fr om that of galaxy 
cluster s. For example, iBaldrv et al.l j2006h and iBolzonella et al.l 
(2009) found a dependence of the mass function on the local den- 
sity. However, as we will show later, comparing field and cluster 
galaxy mass functions, there seem to be no large differences in the 
different environments (see 37. 2\ . 

• star formation: galaxies can increase their stellar mass form- 
ing new stars. The fraction of star-forming galaxies is of course 
much lower in clusters than in the field at the same redshift 
dFinn etall boiO). At least at intermediate redshifts, the star for- 
mation rate in star-forming cluster galaxies is similar to that of 
field star-forming galaxies of the same mass, except for a 10% 
population of cluster star-forming galaxies wi th reduced star for - 
mation rates (Vulcani et al. 2010). For the field, IOeschetalJd2010h 
have found that the latest types of galaxies (bulge-less and irregu- 
lars) show a non-negligible mass growth over the time span from 
z ~ 0.9 to z ~ 0.3, while in earlier types internal star-formation 
since z ~ 0.9 only results in a negligible increase (<1%) in the stel- 
lar mass. 

• Finally, when considering the evolution in the mass functions 
of each morphological type, it is relevant to keep in mind that galax- 
ies can be subjected to a morphological transformation, that is a 
change in their morphological type, either simply due to passive 
evolution or due to influences by other processes. Infalling galax- 
ies meet a new environment that can alter their properties and they 
are subject to several mechanisms that e.g. can quench their star 
formation processes. 

From these descriptions, we can conclude that merging does 
not play an important role in driving the evolution of the mass func- 
tions we observe. 

In principle, harassment could modify the shape of the mass 
distribution, and influence the mass function of each morphologi- 
cal type, altering the galaxy morphology. However, since it is more 
and more efficient at lower masses, we would expect it to decrease 
the number of low mass galaxies at lower redshifts, therefore to 
produce the opposite effect of what we observe (and might have to 
be "overcome" if harassment is an important process). Thus, even 
if it could contribute in shaping and modifying the mass function, 
harassment does not provide a satisfactory explanation for our find- 
ings. 
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Figure 11. Comparison of the mass distribution of EDisCS (blue filled 
squares) and WINGS (black crosses) SO+late-type galaxies. Mass dis- 
tributions are normalized to have the same total number of objects at 
logM»/M0 ^11. Errors are poissonian. Numbers in the labels are the to- 
tal number of galaxies above the respective mass completeness limit. For 
WINGS, both observed (first) and weighted (second) numbers are given. In 
the insets the cumulative distributions are shown. 

The other mechanisms listed above could play an important 
role in driving the observed evolution in clusters, so in the follow- 
ing, we try to quantify separately the extent of the possible contri- 
bution of each one. 



7.1 The role of the morphological transformations 

The infall of galaxies certainly occurs at each redshift, however 
it is hard to quantify its extent. Hence, in this section, we wish 
to test if morphological transformations from one type to another, 
alone, can explain the evolution of the mass function of the various 
morphological types, neglecting any other possible effect, assum- 
ing that the rate of morphological evolution does not depend on 
galaxy mass and considering clusters as closed boxes. In this sec- 
tion we only consider mass-limited samples. 

Since both the mass fraction and the number fraction in ellipti- 
cal galaxies remain quite constant with redshift (as seen also in Ta- 
ble[3j, let us at first assume that the only efficient mechanism is one 
that transforms a late-type into an SO galaxy. If this were the case, 
the mass distribution of late+SO galaxies should not change with 
time. Instead, Figure QT] shows that it is very different at the two 
redshifts (Pks > 99.9%), indicating that the evolution of the mass 
functions cannot be explained simply by turning late-type galaxies 
into SOs and neglecting infall and any other mechanism. 

In the next step, we recreate a more complex situation, includ- 
ing elliptical galaxies in our analysis too and making the hypothesis 
that, still in a closed box, observed early-types in WINGS derive 
from early-types in EDisCS and from a fraction of EDisCS late- 
types that changed their morphology to earlier types. 
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Figure 12. Model of a possible evolution of the mass functions without tak- 
ing into account the infall of galaxies from areas surrounding the clusters. 
Black crosses: observed WINGS early-types, Blue filled squares: observed 
EDisCS early-types (same as third panel Figure [§}■ Red stars: simulated 
mass function of WINGS early-types assuming they derive simply from all 
EDisCS early-types, and from some of the EDisCS late-types (see §6.4 for 
details). 

In practice, since we hypothesize a closed system, we have 
to assume that the fraction of early-types in WINGS (84%) has to 
be equal to the fraction of early-types in EDisCS (47%) plus some 
of the late-type galaxies in EDisCS necessary to maintain the right 
proportions (X times the late-type fraction = X times 58%). We 
determine this quantity assuming that the total number of galaxies 
cannot change and hence X must be equal to 0.64. 

We find (Figure [T2l that the resulting simulated mass func- 
tion is very different from the observed early-type WINGS mass 
function, indicating that the early-types in WINGS cannot derive 
simply from the evolution of late-types of EDisCS, no infall and no 
other mechanism at work. The K-S test support what is clear from 
the figure, giving a probability > 99.99% that the distributions do 
not derive from the same parent one. 

Hence, it is impossible to reconcile the mass functions ob- 
served at low-z with a simple picture in which there is no infall, 
each galaxy type preserves its mass distribution unaltered and late- 
type galaxies are turned into early-types by some mechanism that 
leaves their stellar mass unaltered. 

While morphological transformations are likely to occur, they 
have to be accompanied by other mechanisms to explain the ob- 
served trends. 



In the literature, there are few works in which the mass distri- 
butions of g alaxies of different morphological types are presented. 
Illbert et alj d2010l) analyzed the evolution from z = 2 to z = for 
different morphological and spectral types in the field using the 
COSMOS 2 deg 2 field; in the following, we will consider their 
mass distribution of elliptical and all galaxies in the redshift bin 
0.6 < z < 0.8. Bundy (2005fl analyzed both together and sep- 
arately, for the early- and late-type field galaxies in the GOODS 
fields in the redshift range < z < 1. In the following, we will also 
consider their data between 0.55 ^ z ^ 0.8 as reference. 

First of all, we wish to compare the total mass functions. 
Figure Q~3] shows the comparison between cluster and field mass 
functions for all galaxies, above the common mass limit, that is 
logM»/Mo ~ 10.3 once all masses have been converted to the 
same IMF. The mass functions are normalized in order to have the 
same total number of galaxies above log M* /Mo ^ 1 10 

We see that at high masses (logM*/M0 ^ 11), the mass func- 
tion of field and cluster galaxies at high-z show a rather similar 
shape. However, since we are trying to understand the reasons of 
the observed differences between the mass functions at different 
redshifts in the lower massive tail, where WINGS and EDisCS give 
different results, we are more interested in comparing the shape of 
the mass function of galaxies at intermediate and low masses, as- 
suming the massive tail to be fixed. 

As Figure [T3] illustrates, the two field studies we adopted for 
comp arison give quite different results at these masses. If we fol- 
lowed Ilber t et al.1 ( 120 101) (blue dashed line), we could suggest that 
field galaxies have a steeper mass function than cluster galaxies, 
indicating the presence of a significant environmental mass seg- 
regation. In this way, the infall of galaxies from the field could 
be at least partly responsible for the growth of the low end of the 
mass function with time. Howev er, even if cluster s acquired an in- 
falling population following the Ilber t et all d2O10h mass function, 
this would not fully explain the excess of low mass galaxies ob- 
served in WINGS. The WINGS mass function at low masses is 
too high to be pr oduced by the field mass function observed by 
Illbert et alj d2010l) (see the shape of the WINGS mass function, 
green long dashed line in Figure [TJt. 

Furthermore, Bundy (2005)'s results (red filled hexagons) 
suggest that there are no large differences between the mass dis- 
tribution of galaxies in the different environments at high-z. In this 
case, an environmental mass segregation + infall surely could not 
be responsible for the observed differences between the total mass 
functions in clusters at different redshifts. Unfortunately, based on 
these results, we are not able to securely establish whether field and 
cluster galaxies have similar or different mass distributions. 

Focusing on the different morphological types (plots not 
shown), our preliminary comparison between our cluster sample 
and the field studies suggests that, for all types, the shape of the 
massive end of the mass function does not strongly depend on en- 
vironment. Instead, at lower masses, there is a slight excess of field 
galaxies both for ellipticals and for late-types; while, the compar- 
ison of early-type galaxies (ellipticals + SOs) seems to go in the 
opposite direction, with clusters having proportionally more low- 
mass galaxies than the field. 



7.2 Can infalling galaxies explain the observed evolution? 

Environmental mass segregation in infalling galaxies could play an 
important role, for the reasons discussed in previous sections, but 
we have no means to assess its role without knowing the infall rate, 
mass distribution and morphological mix of accreting galaxies. 



7 From private c om munication, these data are the combination of 
iBundvet alj J2005h and lBundv et~ai] feOOd) 

8 The number of galaxies per unit volume in the field is of course much 
lower than in clusters, but here we are only interested in the shape of the 
distributions. 



18 Vulcani et al. 



oil 



3 



2.5 



1.5 - 



0.5 



XEDisCS 0.4<z<0.8 
-llbert et al. (2010) 0.6<z<0.8 
• Bundy (2005) 0.55<z<0.8 

•WINGS 0.04<z<0.07 
_j i i i I i i i i I i 



10 



10.5 



1 1 
Log, M 



J ♦ 



w 



1 1.5 



12 



0.6 - 



0.4 - 



0.2 




Figure 13. Comparison of the mass distribution of cluster galaxies (EDisCS 
- black crosses) and field galaxi es (llbert & B undy) at similar redshifts . Red 
filled hexagons are from lBundv et alj feOOSl) and lBundv et al . 12006). We 
refer b riefly to them as Bundy (2005); blue dashed lines are from llbert et al. 
120101) . For comparison, we also plot the Schechter fit of the WINGS mass 
function (green long dashed line). Mass distributions are normalized to have 
the same total number of objects at logM* /Mq ^11, Errors are poissonian. 



We have to note that these are only preliminary results, per- 
formed using heterogenous data and slightly different redshift 
ranges, so it cannot be used to draw definite conclusions. In the 
future, to establish the existence and quantify the eventual extent of 
an environmental mass segregation, it will be important to compare 
mass functions of different environments at the same redshift using 
homogenous data. To understand the evolution of the mass func- 
tion in clusters, it will be especially useful to consider the mass 
function in the cluster outskirts, whose galaxies will have time to 
become part of clusters before z = 0. 

Nevertheless, this preliminary comparison between the cluster 
and field mass functions at high-z seems to suggest that the cause 
for the evolution of the mass function in clusters probably needs 
to be sought elsewhere than in environmental mass segregation ef- 
fects. 



7.3 How important is star formation? 

We can attempt to evaluate the importance of the contribution of 
star formation in changing the mass function with redshift, know- 
ing the fraction of star-forming galaxies and the rate of mass growth 
due to star formation. This process doesn't act directly by changing 
the morphological type, but it increases the galaxy st ellar mass. As 
previously demonstrated by several works (see e.g. Noeske et al.l 
2007; Feulner et al. 2005; Perez-Gonzalez et al. 2005 and our own 



Vulcani et al.ll20irj results), low-mass galaxies have higher values 



of specific star formation rate (hereafter SSFR, star formation rate 
per unit of galaxy stellar mass) than high mass galaxies so they 
proportionally increase their mass much more rapidly than more 



Figure 14. Fraction of blue galaxies over all galaxies as a function of the 
mass in the EDisCS mass limited sample. Galaxies with (U- V) rea _f rame ^ 
0.8 are considered blue, hence star-forming. The red line represents the least 
squares fit performed on the points. 



massive galaxies. As a consequence, there can be some galaxies 
whose mass is below the mass limit at the redshift of EDisCS, but 
with time can grow so that they enter the sample at the redshift of 
WINGS, changing the shape of the mass function at low- interme- 
diate masses. 

We build a toy model to try to quantify the importance of star 
formation. First, we discuss an extreme case. We "simulate" a pos- 
sible situation in which infall does not act and the only event is 
the increase of mass through star formation occurring in galaxies 
already located in clusters. 

In clusters, star formation cannot act undisturbed from z ~ 0.7 
to the current epoch, since there are probably several processes that 
cause star formation to end. Therefore, we simulate a situation in 
which star formation is left undisturbed only for 3 Gyr and then it 
stops suddenly. During the 3 Gyr, we hypothesize that it is not in- 
fluenced by some external processes, so it decreases naturally with 
time, according to the evolution with redshift of the SSFR-Mass 
relation. We have chosen 3 Gyr as a "generous" timescale, as this 
can be considered a sort of upper limit for the time a galaxy will 
continue forming stars undisturbed in clusters, given the dynamical 
timescale in clusters (~ 1 Gyr) and the long timescale expected for 
some of the quenching processes proposed in clusters (e.g. stran- 
gulation, a couple of Gyrs)0 

First of all, we take into account the fact that only star-forming 
galaxies grow in mass due to star formation, while passive galax- 
ies do not. We consider as star-forming only EDisCS galaxies with 
(U — V)rest-frame ^ 0.8 (blue galaxies). For each mass, we estimate 
the fraction of star forming galaxies, calculated as the number of 

9 Even allowing more time for the SFR to act (i.e. 6 Gyr), we found that 
star formation in galaxies that are already in clusters at high-z cannot be the 
only mechanism responsible for the observed evolution (see below). 
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blue galaxies over all galaxies (Figure [14}. Since the relation be- 
tween the fraction of blue galaxies and logM seems to be linear, 
we perform a least squares fit and we extrapolate the fraction also 
at masses below the mass completeness limit. When the fraction 
would t urn out to be greater than 1, it is just set equal to 1. 

In lVulcani et al.U2oToh we analyzed the trend of the SSFR in 
clusters as a function of galaxy mass. We found that this trend is 
very similar to the trend found by Noe ske et alj {2007) for the field. 
The main difference is that in clusters there is a population of star 
forming galaxies (~ 10%) that has a reduced star formation rate for 
their mass and that probably will e nd the process of sta r formation 
in the near future (for details, see I Vulcani et So, when 

we consider the fraction of star forming galaxies we eliminate from 
the star-forming fraction 10% that corresponds to galaxies with re- 
duced star formation. 

Since we need the value of median SSFR at low masses and 
our cluster mass completeness limit is quite high, to estimate the 
mass gr owth we use the med ian SSFR-Mass relation found for the 
field by iNoeske et al] d2007h . calculated in mass bin of 0.1 dex. 
To reach even smaller masses, we extrapolate this relation down to 
logM./Mo ~ 8.1 (the minimum mass of a galaxy at z ~ 0.7 that 
it turns out at z ~ could enter our sample due to star formation). 
We estimate by how much galaxies can increase their mass in 3 
Gyr from z ~ 0.7 adopting the SSFR appropriate for the epoch and 
the mass considered. In our estimates, we account for the fact that 
the observed stellar mass at any time is equal to the integral of the 
mass of stars formed before, minus the fraction of mass returned 
by stars into the interstellar phase (about 25-35% of the total mass 
every time the galaxy doubles its mass). 

Finally, we assume that the EDisCS total mass function at 
masses below our mass limit remains flat, i.e. has a similar behavior 
to what we observe between logM, /Mq = 10.2 and logM»/M0 = 
10.8. To verify this assumption where we could, we considered blue 
star-forming galaxies, for which we can reach a mass completeness 
limit lower than that of all galaxies, that is logM* /Mq = 10.0. Us- 
ing the relation given in Figure[l4] we then estimate the total num- 
ber of galaxies at these low masses and we find the result to be 
consistent with our hypothesis that the EDisCS total mass function 
continues flat to low masses, at least down to logM„/Mo = 10.0. 

According to the fractions of star forming galaxies, we add 
in each mass bin of the observed EDisCS total mass function, the 
number of galaxies that could enter that bin at z ~ 0, and we sub- 
tract the number of galaxies that left that bin with time. Then we 
normalize the new mass function so that the number of galaxies at 
high masses (logM,/MQ = 11.5, where the growth turns out to be 
negligible), is the same as that in the EDisCS mass function. 

Figure Q3] shows the results of our attempt: we find that high- 
z cluster members whose mass is increased due to star formation 
could only partially fill the observed gap between WINGS and 
EDisCS at low/intermediate masses. In fact, this test could fully 
account for the growth of the mass function only in the first bins, 
for logM*/MQ 10.4, while the growth of intermediate masses 
still remains unexplained. 

We stress that this result is found assuming that the high-z 
mass function is flat at masses below our limit, that infall is negli- 
gible and that star formation continues undisturbed for only 3 Gyr 
in the cluster environment. 

Then, as a second more realistic test, we also take into ac- 
count the contribu tion of infalling galaxies from the field, using 
lBundvetalJ ( l2005h data. Since their data are given in mass bins of 
0.3 dex, we use their same binning. As we do not know the infall 
rate, we cannot estimate the exact relative contribution of cluster 
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Figure 15. Comparison of the mass distribution of WINGS (black crosses) 
and EDisCS modified under the assumption that star formation produces 
an increase at intermediate masses ("EDisCS grown", red empty stars). For 
reference, also the EDisCS mass function is plotted (blue filled squares). 
WINGS and EDisCS mass functions are normalized as usual. The mass 
distribution of EDisCS grown is normalized to have the same total number 
of objects at logM, /Mq > 11.5 as EDisCS. Errors are poissonian. Num- 
bers in the labels are the total number of galaxies above the respective mass 
completeness limit. For WINGS, both observed (first) and weighted (sec- 
ond) numbers are given. 



and field galaxies. So, we normalize both the total mass function 
of EDisCS and that of the field to 1 and we sum them together, 
giving a double weight to the field galaxies. This correspond to the 
amount of cluster mass growth between the two redshifts from dark 
matter simulations (Table 4 in Poggianti et al. 2006). Then we as- 
sume that the star forming fraction at each mass is the same as that 
used in the previous test[^j To extrapolate the field mass function 
at lower masses than the Bundy's completeness limit, we perform 
a Schechter fit of the field mass function. To do this, we use only 
Bundy's late-type galaxies since we assume that the great majority 
of star-forming galaxies belongs to this morphological type. After- 
wards, in agreement with the previous toy model, we compute the 
mass growth of galaxies (both of those already in clusters and of 
those infallen) due to star formation in the lapse of time of 3 Gyr 
and estimate the new "grown" mass function. 

Results are shown in Figure [T6] Comparing the shape of the 
"grown" mass function, we find that it is more similar to the 
WINGS mass function than to the EDisCS one, from which it is 
derived. Our result indicates that the SFR can indeed give a non- 
negligible contribution to the evolution of the mass function. In 
this test, infalling galaxies are able to fill the gap also at interme- 
diate masses, while this did not occur in the previous test. Impor- 



10 We can't estimate the real star forming fraction - mass relation in the 
field, hence we consider the relation we have found for clusters. Certainly, 
this is a hard lower limit and allows us to be as conservative as possible. 
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Figure 16. Comparison of the mass distribution of WINGS (black crosses) 
and EDisCS modified considering the contribution of the infalling galax- 
ies and under the assumption that SF increases the mass of intermediate 
mass galaxies ("EDisCS grown", red empty stars). For reference, also the 
EDisCS mass function is reported (blue filled squares). Mass functions are 
normalized as usual. Errors are poissonian. 



tantly, the shape of the new mass function is more similar to that of 
WINGS than to EDisCS. 

We emphasize that in our model we have used for the field 
the star-forming fraction versus mass relation observed for clusters. 
This is surely a conservative lower limit, and therefore results in a 
lower limit for the contribution of infallen galaxies. We also point 
out that we are supposing that all field galaxies enter the clusters at 
z ~0.7, while galaxies enter at different epochs, and those entering 
later will have more time to continue forming stars and contribute 
to the growth of the mass function. 

Moreover, we do not, for example, consider the effect of the 
enhanced star formation (starb ursts) that might to be associated 
with galaxy infall into clusters jDressler et alj[20()j^) . This would 
provide an additional mass growth that helps reconcile the EDisCS 
and WINGS mass functions. 

From these models, we suggest that star formation in both 
cluster galaxies and, most of all, in later infalling galaxies is proba- 
bly the dominant mechanism responsible for the evolution of the to- 
tal mass function. We have seen that star formation in star-forming 
galaxies that are already in clusters at high-z does not fully explain 
the observed gap in the mass functions at the two redshifts, while 
adding the infalling field population more closely reproduces the 
observed mass function at low-z. 

This shows that infall does play a role, contributing a large 
number of low mass galaxies that, as long as they continue form- 
ing stars, increase their mass due to star formation. In §7.2 we have 
seen that it is currently unclear whether, in addition to this effect, 
infall also contributes to the evolution of the mass function via en- 
vironmental mass segregation, by directly increasing the relative 



numbers of low-mass (logM»/MQ = 10.2- 10.8) versus high-mass 
galaxies logM»/M ^ 10.8. 

In our test we ha ve conservatively adopted the mass function 
of lBundvet"al] j2005) that, as we have seen in §7.2, is more similar 
to the EDisCS mass function and suggests no environmental mass 
seg regation. If t he fiel d mass function was steeper at lower masses 
(as lllbert et ail {2010) find), the contribution of the mass growth 
due to star formation in infalling galaxies could be even stronger. 

Obviously, star formation can affect the evolution of the late- 
type galaxy and total mass functions, but alone cannot explain 
the observed evolution of the mass function of early-type galax- 
ies. Moreover, star formation cannot be directly responsible for the 
observed evolution in the morphological mix. Some other mecha- 
nism/s must be at work transforming late-type galaxies into early- 
types, thereby changing the shape of the early-type mass func- 
tionsE] 



8 SUMMARY AND CONCLUSIONS 

In this paper, we have characterized the stellar mass functions and 
their evolution from z = 0.8 to z = of two mass limited samples 
of galaxies in clusters drawn from the WINGS (logM„/Mo > 9.8) 
and EDisCS (logM*/M > 10.2) datasets. We considered both the 
total mass function and the mass function of different morpholog- 
ical types. We then analyzed the morphological fractions both for 
mass limited and magnitude limited samples. 
Our main results are: 

• The total stellar mass distribution of galaxies in clusters 
evolves with redshift. In a mass-limited sample, clusters at high- 
z have on average more massive galaxies than at low-z. Assum- 
ing no evolution at the massive end, the population of 10.2 ^ 
log M* /Mo ^ 10.8 galaxies must have grown significantly between 
z = 0.8 andz = 0. 

• At low redshift, different morphological types have signifi- 
cantly different mass functions, although SOs and late-type galaxies 
more massive than logM»/MQ = 10.2 may have similar mass dis- 
tributions. At high redshift our results are not conclusive. 

• The mass distribution of each morphological type evolves 
with redshift. In all types there are proportionally more massive 
galaxies at high- than at low-z. 

• Both for the mass-limited and the magnitude-limited sam- 
ples, the proportion of galaxies of different morphological types as 
a function of galaxy stellar mass (the morphology-mass relation) 
evolves strongly with redshift, as expected given the evolution of 
the morphological fractions. At both redshifts, in the mass limited 
sample, the incidence of each type does not strongly depend on 
mass up to log M, /Me ~ 11. Above this value, ellipticals domi- 
nate. 

• Both considering a mass limited and a magnitude limited 
sam ple, and in agreement with previous magnitude-limited works 
(e.g. iDressler et al.ll 19971; iFasano et alj|200ol: iPostman et alj|2005l: 
ISmith et al.ll2005l : lDesai et alj2007MPoggianti et al.ll2009h . we find 
that the number fraction of elliptical galaxies is almost constant 



1 1 We note that such mechanism can in principle be itself mass-dependent, 
and could transform low-mass galaxies more efficiently than high-mass 
galaxies. In this scenario, also harassment, by stripping mass to galaxies 
and making them change their morphology, could partly contribute to the 
change of the mass functions type by type. 
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with redshift (~ 30%). In contrast, the fraction of SOs and late- 
types considerably changes with time (SOs ~ 50% at the current 
epoch, late-type ~ 55% in the distant Universe), as also reported in 
those other studies. 

• In addition to the number fraction of ellipticals, also the frac- 
tion of stellar mass in ellipticals remains constant with redshift, 
being about 40% both in EDisCS and in WINGS. At high-z, 13% 
of the stellar mass is in SOs, and 46% in later types. At low-z, the 
mass fractions are inverted: 42% in SOs, and 15% in late-types. 

The most likely explanation of the observed evolution of the 
mass functions is the mass growth of galaxies due to star formation 
in both cluster galaxies and, most of all, in galaxies infalling from 
the cluster surrounding areas. In this way, galaxies that are not part 
of our distant universe sample, after a reasonable lapse of time, can 
increase their mass due to star formation so that they can enter our 
local Universe sample. 

Therefore, infall does play a role contributing a large num- 
ber of star-forming galaxies growing their mass. In principle, infall 
could also contribute by directly increasing the population of low- 
mass galaxies compared to the massive population, due to infalling 
galaxies having a different mass distribution from galaxies already 
in clusters at high-z. However, we do not find conclusive evidence 
for such environmental mass segregation in our comparison with 
field studies. 

In the future, it will be important to compare the galaxy stellar 
mass distribution in cluster cores, cluster infalling regions, groups 
and the field using homogeneous data, to conclusively establish 
whether and how much the total galaxy stellar mass function de- 
pends on environment at different redshifts. 

Finally, star formation, acting mainly on late-type galaxies, 
can directly affect the evolution of the total and late-type galaxy 
mass functions, however it cannot directly account for the observed 
evolution of early-type galaxies. Subsequent morphological trans- 
formations from late to early types likely cause the evolution of the 
mass functions of SOs and ellipticals. 
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APPENDIX A: THE RELIABILITY OF EDISCS PHOTO-Z 
MASS FUNCTIONS 

Throughout the paper, for EDisCS, we have used the photo-z mem- 
bership. Here we show that, in the mass range in common, the mass 
function determined from photo-z's is in agreement within the er- 
rors with the mass function determined using only spectroscopic 
members and spectroscopic completeness weights, derived simi- 
larly to what we have done for WINGS. The two mass functions for 




Log, M 

Figure Al. Comparison between the mass distribution of all EDisCS galax- 
ies with spectroscopic membership (black crosses) and photo-z membership 
(blue filled squares). Mass distributions are normalized to have the same to- 
tal number of objects at logM»/A/Q > 11. Errors are poissonian. In the 
insets the cumulative distributions above the same mass limit are shown. 

the mass limited sample are shown in (Figure lATt . In this compar- 
ison, we are using all photo-z and spectroscopic members, without 
considering their i?200 or their morphology. 

A K-S test above \ogM* / Mq = 10.6 cannot reject the hypoth- 
esis that the two samples are drawn from the same parent distribu- 
tion (K-S probability ~ 82.1%). This gives us confidence in using 
the photo-z's, that allow us to reach much deeper mass limits. 



APPENDIX B: THE MAGNITUDE COMPLETENESS 
LIMIT 

In this Appendix we show the mass distributions in magnitude lim- 
ited samples, considering the same magnitude limits adopted by the 
morphological fraction studies (§6.1). We show how this choice al- 
ters the final results: in fact, these mass distributions are very dif- 
ferent from those we derive using mass limited samples. 

In mass selected samples only a modest and systematic evolu- 
tion of galaxy masses with time is expected, unlike magnitude se- 
lected samples where galaxies can enter and leave the sample due 
to variations in their star formation activity. A mass selection gives 
a more robust mean of identifying likely progenitors of z=0 cluster 
galaxies in higher redshift clusters. However, we note that also in 
mass limited samples some progenitors could lack. First, galaxies 
can continue to grow in mass and so some galaxies that weren't in- 
cluded at high redshift enter in the sample at low redshift, moreover 
the infall of galaxies from the field can alter the mass distribution. 

It is important to note that the choice of a magnitude limit 
implies a natural mass limit below which the sample is incomplete. 
This corresponds to the mass of the reddest galaxy at the limiting 
magnitude. 
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all 




9.5 10 10.5 11 11.5 

Log, M 



Figure Bl. Comparison of the mass distribution of EDisCS (blue filled 
squares) and WINGS (black crosses) galaxies, of all morphological types, 
for the magnitude limited samples. Mass distributions are normalized to 
have the same total number of objects at log M* / Mq ^11. Errors are pois- 
sonian. In the insets the cumulative distributions are shown. 



For WINGS, M v = -19.5, and B- V = 1.2, the mass complete- 
ness limit corresponds toaM~ 10 10 ' 5 M@. 

For EDisCS, M v = -20, and B- V = 0.9, it corresponds to a 
M~ lO la4 M . 

When using magnitude limited samples, it should therefore 
always be kept in mind that below the mass completeness limit the 
sample is incomplete, hence results will differ from those obtained 
from a complete (e.g. mass limited) sample. In the following, we 
show some examples. These plots have to be directly compared 
with those already presented in this paper. 

Figure IbTI shows the mass distribution of all galaxies, regard- 
less of their morphological type, for WINGS and EDisCS galaxies 
in the magnitude limited sample. Both distributions rise, reach a 
peak, then decline again. 

For WINGS, for which the completeness mass limit of the 
magnitude limited sample is much higher than the limit of the 
mass limited sample, it is clear that the shape of the distribution 
below the completeness mass limit is very different from the one 
obtained for the mass limited sample (cfr. Figure IbTI with Fig. 3). 
Instead, no significant change is expected nor observed between 
the mass limited and the magnitude limited distributions above 
logM„/M Q ~ 10.5. 

As for the mass limited samples, the mass functions of 
WINGS and EDisCS are very different. If we consider the whole 
mass range covered by galaxies in the magnitude limited sample, 
a K-S test rejects the null hypothesis of similar mass distribution 
with a probability 100%. 

Then, we also show (Figure lB2t the comparison of the mass 
distribution of the different morphological types at the two red- 
shifts, to be contrasted with the mass-limited ones in Fig. 5. Again, 
in the magnitude limited samples, incompleteness setting in below 



the completeness mass limit creates a spurious declining trend at 
low masses. 

The comparison of magnitude-limited and mass-limited sam- 
ples clearly shows that the mass distribution derived from the for- 
mer is meaningless, because affected by incompleteness, below the 
the limit corresponding to the mass of a galaxy with the reddest 
color and the faintest magnitude in the sample. 
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Figure B2. Comparison of the mass distribution of EDisCS (blue open circles) and WINGS (red stars) galaxies, for different morphological types, in magnitude 
limited samples. Mass distributions are normalized to the total number of objects with logM»/MQ ^ 11. In the insets the cumulative distributions are shown. 
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